hep-ph/9905403 



CERN-TH/99-109 
DESY 98-199 



o> 

On 
On 

3 



ON 



Measuring Gauge- Mediated SuperSymmetry Breaking 
Parameters at a 500 GeV e + e Linear Collider* 

Sandro Ambrosanio a c and Grahame A. Blair 6 c 

a CERN - Theory Division, 
CH-1211 Geneva 23, Switzerland (current address) 
e-mail: ambros@mail.cern.ch 



b Royal Holloway and Bedford New College, 
University of London, Egham Hill, Egham, 
Surrey TW20 OEX, U.K. (current address) 



e-mail: g.blair@rhbnc.ac.uk 



c Deutsches Elektronen- Synchrotron DESY, 
;> NotkestraBe 85, D-22603 Hamburg, Germany 

m ' 
o 
^t- 
m 
o 

ON • Abstract 
ON 

\ We consider the phenomenology of a class of gauge-mediated supersymmetry (SUSY) breaking 

(GMSB) models at a e + e~ Linear Collider (LC) with -E c . .m. up to 500 GeV. In particular, 
we refer to a high-luminosity (£ ~ 3 x 10 34 cm -2 s" 1 ) machine, and use detailed simulation 
tools for a proposed detector. Among the GMSB-model building options, we define a simple 
framework and outline its predictions at the LC, under the assumption that no SUSY signal 
is detected at LEP or Tevatron. We assess the potential of the LC to distinguish between the 
various SUSY model options and to measure the underlying parameters with high precision, 
including for those scenarios where a clear SUSY signal would have already been detected at 
the LHC before starting the LC operations. Our focus is on the case where a neutralino (Ni) is 
the next-to-lightest SUSY particle (NLSP), for which we determine the relevant regions of the 
GMSB parameter space. Many observables are calculated and discussed, including production 
cross sections, NLSP decay widths, branching ratios and distributions, for dominant and rare 
channels. We sketch how to extract the messenger and electroweak scale model parameters from 
a spectrum measured via, e.g. threshold-scanning techniques. Several experimental methods 
to measure the NLSP mass and lifetime are proposed and simulated in detail. We show that 
these methods can cover most of the lifetime range allowed by perturbativity requirements and 
suggested by cosmology in GMSB models. Also, they are relevant for any general low-energy 
SUSY breaking scenario. Values of ct^ as short as 10's of um and as long as 10's of m can 
be measured with errors at the level of 10% or better after one year of LC running with high 
luminosity. We discuss how to determine a narrow range ( < 5%) for the fundamental SUSY 
breaking scale y/~F, based on the measured m ^ , CTjy . Finally, we suggest how to optimise the 
LC detector performance for this purpose. 
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1 Introduction 



If the world is supersymmetric at short distances, then the gauge hierarchy problem can 
be naturally solved. The most compelling proof of this hypothesis would be direct de- 
tection of superpartners at colliders. This has not been achieved so far, which tells us 
that supersymmetry (SUSY) must be broken. In order for a SUSY theory to preserve 
its theoretically pleasant characteristics, supersymmetry breaking (SSB) can only occur 
in a "soft" way []T|. However, this constraint still allows a general phenomenological ap- 
proach to SSB involving over a hundred new parameters in addition to the Standard 
Model (SM) ones. Strategies for searches at present and future colliders must then rely, 
at least to start with, on theoretically well-motivated schemes for SSB, providing a more 
definite framework and living on a manageable parameter space. A related question is 
how this SSB is transmitted to the visible (light) sector of the theory, e.g. the particles 
of the Minimal SUSY extension of the SM (MSSM). Historically, the most popular ap- 
proach has been that SUSY is broken at very high energies (HESB) of the order of the 
Planck mass or the scale of Grand-Unified Theories (GUT) and SSB is communicated 
to the MSSM sector through gravitational interactions. Such an approach goes usually 
under the name of (minimal) Supergravity [(m)SUGRA] or, with some additional as- 
sumptions, Constrained MSSM (CMSSM) H. More recently, another equally attractive 
scenario has earned large consensus and recognition, both among theorists and experi- 
mentalists, the Low-Energy Supersymmetry Breaking (LESB) option, and in particular, 
the Gauge- Mediated (GMSB) version of it [Q]. LESB, in itself, may already have striking 
phenomenological consequences, as it was shown in pioneering works by Fayet j3J. Indeed, 
gravity enters the expression for the gravitino mass, 

m 3/2 =mf= ^ ~ ( ^ \ 2.37 eV, (1) 

where yF is the fundamental scale of SSB, 100 TeV is a typical value for it in LESB 
models, and M' P = 2.44 x 10 18 GeV is the reduced Planck mass. As a result, the gravitino 
is so light in LESB models that it plays always the role of the lightest SUSY particle 
(LSP) and can be treated as massless for all kinematics purposes at high energy colliders. 
However, for \f~F~ <C M' P) the dominant gravitino interactions come from its longitudinal, 
spin-1/2 components, namely the goldstino components that the gravitino has acquired 
through the so-called SUSY-Higgs mechanism. Hence, gravity does not enter the strength 
of the gravitino couplings to matter, which in the relevant approximation are proportional 
to the mass splitting between superpartner masses and the ordinary SM particle masses 
and inversely proportional to F. The latter can be small enough to render the gravitino 
relevant for collider phenomenology. (It should be noted here that a light gravitino LSP 
can also be obtained within the framework of no-scale SUGRA models ||.) 

The phenomenological scenario in LESB with conserved -R-parity (which we assume 
in the rest of the paper) can be summarised as follows: 

• every produced SUSY particle has to decay to the G, possibly through a cascade; 

• since the goldstino interactions are still much weaker than the ordinary SM gauge 
and Yukawa interactions, every decay chain has to involve the next-to-lightest SUSY 
particle (NLSP), which in turn will finally decay to the gravitino; 

• depending on \fF, the production energy and details of the SUSY spectrum, the 
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NLSP can decay close to the interaction point (i.p.), within or outside a collider 
detector, producing a plethora of new spectacular signatures. 

Among the possible mechanisms for transmitting LESB to the MSSM fields, by far the 
most effective and theoretically satisfying is GMSB, where a so-called messenger sector 
is responsible for communication between the secluded sector where SSB takes place and 
the visible sector, via SM gauge interactions. Mainly motivated by a natural suppression 
of the SUSY contributions to flavour- changing neutral current (FCNC) and CP-violating 
processes, such a scenario was first explored in various forms in several early 1980 's works 
|J and then recently revived in its present version in the famous papers of Ref. jTJ. 
Remarkably, in addition to the appealing theoretical features, the minimal version of 
GMSB also provides a powerful tool for building very predictive models and calculating 



spectra from just a handful of parameters, as done e.g. in Refs. M,M, 11]. An important 



boost to the popularity of LESB and GMSB models came a few years ago due to a possible 



explanation of the anomalous CDF e + e~ / y / yE/r event within this framework ]12,[13|. Today, 
such an explanation seems more unlikely, yet it worked fine in stimulating a considerable 
number of dedicated analyses and searches for GMSB-inspired new signals at LEP and 
Tevatron, which are of course of much broader interest |12|, [14], [L5|. Hence, it is now time 
to think about how similar searches could be pursued at next generation colliders and 
how the reach in the GMSB parameter space of such machines could be optimised. Some 
work in this respect has already been carried out for the Tevatron Run II jn| and for the 



LHC [ lTf . In this paper, we will be instead mainly concerned with GMSB phenomenology 



at a first phase of operations of a e + e~ Linear Collider (LC) with c.o.m. energy up to 
around 500 GeV, and will focus on the case where a neutralino is the NLSP. In particular, 
we will refer to a high-luminosity machine with C ~ 3 x 10 34 cm -2 s _1 , such as being 
considered e.g. by the ECFA/DESY TESLA project, and the related proposed detector. 
Many of our results and experimental methods can be easily extended to more general 
LESB models and might even have an impact on other scenarios such as HESB models 
with .R-parity violation, where delayed NLSP or LSP decays can take place. 

The rest of this paper is organised as follows. In Sec. 0, we briefly describe our GMSB- 
model building framework and specify the region of the parameter space we are interested 
in here. In Sec. ^ we focus on the general phenomenology of models with a neutralino 
NLSP and discuss its possible (delayed) decays, including some new aspects of interest for 
the LC. In Sec. (|, we introduce the main features of the proposed TESLA linear collider 
and give the expected machine parameters relevant to our study. In Sec. ||, we discuss 
the general characteristics of the GMSB signal at the LC and show an example of how 
it is possible to extract a good amount of information about the GMSB parameters via 
a simple experimental technique in principle possible at such a machine. In Sec. |6], we 
describe the relevant characteristics of the LC detector and the software we used for our 
simulations. In Sec. [7], we discuss several methods for measuring the neutralino NLSP 
properties, and in particular its mass and lifetime, using different parts of the detector, 
and we show our results. Finally, in Sec. |8|, we draw our conclusions and comment on 
how the performance of the LC in measuring GMSB parameters depends on details of the 
machine and detector design. We also give a few suggestions to optimise such performance. 



2 Models with Gauge- Mediated SUSY Breaking 

In addition to the automatic suppression of SUSY FCNC, GMSB models have many 
other interesting characteristics. For instance, the sparticles' masses have a transparent 
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and common origin and all approximately scale with a single parameter A which is the 
universal soft SUSY breaking scale for the visible sector. Also, the resulting spectrum is 
notably different from other SUSY scenarios. Further, it is possible to achieve radiative 
electroweak symmetry breaking (EWSB) nicely. There are however problems connected 
for instance with the lack of a compelling dynamical mechanism for generating the SUSY 
parameter /z, but this is common to other SUSY frameworks. 

As far as GMSB-model building is concerned, we will follow closely the approach used 



in Ref. flll[] for LEP2 phenomenology, with some extensions of the parameter space to 
account for the wider kinematical reach of a LC. We will not repeat the technical details 
here, but in order to fix our framework and notations we remind that, after imposing 
EWSB, a minimal GMSB model can be constructed from the following parameters, 

M mess , A" mess , A, tan/3, sign(/x), (2) 

where M mess is the overall messenger scale; N mess is the so-called messenger index that 
parameterises the structure of the messenger sector; A is the universal soft SUSY breaking 
scale felt by the low-energy sector; tan/3 is the ratio of the vacuum expectation values 
(VEVs) of the two Higgs doublets; sign(/i) is the ambiguity left for the SUSY higgsino 
mass after EWSB conditions are imposed. The MSSM parameters and the sparticle 
spectrum are determined from renormalisation group equation evolution starting from 
boundary conditions at the M mess scale, where M a = N mess Ag(A/M mess )a a , (a = 1, 2, 
3) for the gaugino masses and rh 2 = 2A" mess A 2 /(A/M mess ) J2a{ a a/4=7i) 2 C a for the scalar 
masses. Here g, f are the one, two-loop functions whose exact expression can be found 
e.g. in Ref. fllR , and C a are the quadratic Casimir invariants for the scalar fields. The 
Af couplings are taken to be zero at the messenger scale, since they are generated (first 
power) at the two-loop level. We use a phenomenological approach for Bfj,, which is not 
assumed to vanish at M mess , but is instead determined together with by requiring 
correct EWSB. 

For the purpose of exploring the GMSB parameter space of interest for the LC, we 
generated about 20,000 models, of which about 5,000 have a neutralino NLSP. The spec- 
tacular GMSB signatures, most of which are free from SM-background, make it generally 
possible to exclude GMSB models at LEP2 with ttinlsp < few GeV We esti- 



mate that in a few years searches at LEP and Tevatron will only allow models where the 
whole MSSM spectrum is above about 100 GeV, at least in most typical GMSB scenarios. 
(A remarkable exception is the case where the neutralino is the NLSP and decays outside 
the detector, due to relatively large values of a/F, but this is of no special interest here.) 
Hence, we limit ourselves to models where 100 GeV < ttinlsp < 250 GeV = t/s hC /2. As 
a result, the relevant range for A is between about 60 TeV/A" m e SS and 200 TeV / \f N mess . 

For the sake of simplicity, at first we considered only models where N mess is a positive 
integer between 1 and 10 (actually, we could not construct a model with A" mcss > 8 
satisfying all constraints described above and below). As an example, if the messenger 
sector consists of a 5+5 of the global GUT group SU(5) D SU(3)c®SU(2) L ®U(l)y, then 
A^mess = 1, while if it also includes a 10+10, then A^ss = 1 + 3 = 4. Our messenger scale 
M mess is bounded from below by several constraints. First, to avoid excessive fine-tuning 
of the messenger masses, we impose M mess > 1.01A. Second, we require that the mass 
of the lightest messenger scalar be much heavier than the MSSM particles (at least 10 

TeV, that is M mess > A +V A2 +( 20 TeV ) 2 ,). Finally) to 

preserve gauge-coupling unification, 
we also impose M mess > Mqut exp(— 125/A" m e SS ). In this way, the lowest allowed value 
we obtained for M mess is around 19 TeV. Further, to start with, we set a nominal upper 
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bound on the messenger scale M mess < 10 5 A =>■ M mcss < 2 x 10 10 GeV. We will see 
that this is overruled by other constraints described below. As for tan (3, we require it to 
be larger than 1.2 (to avoid imminent bounds from SUSY Higgs searches at LEP2 and 
non-perturbative blowing up of the top Yukawa coupling below the GUT scale) and we 
could not construct a coherent model with correct EWSB and tan (3 larger than about 55, 
with a mild dependence on A. 

In addition to these parameters, for each given GMSB model, a value for the funda- 
mental SUSY breaking scale y/~F has to be specified to complete the information needed 
for collider phenomenology. The ratio F s /F = AM mess /F, where F s is the scale of SUSY 
breaking felt by the messenger particles, depends on details of the secluded sector and the 
communication between it and the messengers. If this occurs, e.g., via a direct interac- 
tion and the goldstino superfield coincides with a single superfield entering the messenger 
superpotential (which we will assume in the following for simplicity), then one can infer 
from perturbativity arguments up to the GUT scale that the corresponding coupling has 



to be smaller than one ]TT|]. In models with radiative secluded-messenger communication, 



the ratio can be even much smaller. In general, one can argue that 



\J~F > \J AM mcss > A. (3) 

This allows the determination, for each given GMSB model, of a lower bound for the 
gravitino mass (and the NLSP lifetime, as we will see) and an upper bound for the 
strength of its interactions with matter ~ 1/F. In our set of models of interest for the 
LC with 100 GeV < m NLS p < 250 GeV, we find m & > 0.2 eV and VF > 30 TeV. 
Unfortunately, there is no such compelling argument to put a strict upper limit on \J~F~ 
that can be of relevance to collider physics. In a simple cosmological scenario, one might 
invoke the argument that if the gravitino mass is too heavy ( > 1 keV ==>- \[F > 
few thousand TeV), then the gravitino relic density could over-close the universe fl8| . 



This is of some use for our purposes and we will exploit this argument in the following. 
However, one has to keep in mind that a heavier gravitino can well be in agreement with 
cosmological scenarios including an inflationary epoch. Barring the latter possibility, one 
finds that the upper limit on the gravitino mass can only be satisfied in our framework if 
-^mess ^ 2 x 10 8 TeV, in models of interest for the LC. This also implies that values of 
A^mess larger than 6 are highly disfavoured in this case. 

In this parameter space, we generated models by means of a private computer program 
called SUSYFIRE [B5|, an updated, generalised and Fortran-linked version of the program 



used in Ref. ||11|| , which can produce minimal and non-minimal GMSB and SUGRA 
models. For scanning, we used logarithmic steps for A, A/M mess and tan/3. The program 
proceeds by iterating the following: setting the masses and the gauge couplings at the weak 
scale; evolving the RGE's to the messenger scale; setting the messenger scale boundary 



conditions (see Eqs. (23), (24) in Ref. [IT]) for the soft sparticle masses; evolving the RGE's 
back to the weak scale, taking care of decoupling each sparticle at the proper threshold. 
We use two-loop RGE's for the gauge couplings, third generation Yukawa couplings and 
gaugino soft masses. The other RGE's are at the one-loop level. We require EWSB using 
the one-loop effective potential approach (one-loop Higgs masses + consistent corrections 
from stops, sbottoms and staus) at the ^m^m^ scale and we eliminate and Bfi in 
favour of tan f3 and Mz- 

The phenomenology of GMSB models is largely dependent on which particle is the 
NLSP or, better, on which sparticle(s) has (have) a large branching ratio (BR) for decaying 
to its SM partner and a gravitino. Four main scenarios are possible: 
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Neutralino NLSP scenario: Occurs whenever < (m^ — m T ). Here typically a 
decay of the Ni to G7 is the final step of decay chains following any SUSY production 
process. As a consequence, the main inclusive signature at colliders is prompt or 
displaced photon pairs + X + missing energy Ni decays to GZ° and other minor 
channels are also important for this study, as we will see in the following. In the 
rest of this paper, we will focus on this possibility, although we are well aware that 
the other scenarios are very relevant for LC phenomenology and we plan to devote 
further work to them. A detailed discussion of the neutralino NLSP case will be 
carried out in Sec. [| 

Stau NLSP scenario: Defined by m? 1 < Minfm^, Trig ] —m T , features f\ — > Gr decays, 
producing r pairs or charged semi-stable f\ tracks or decay kinks + X + missing 
energy. Here £ stands for e or [l. 

Slepton co-NLSP scenario: When mg < Minfm^, + m T ], £r — > Gi decays are 
also open with large BR. In addition to the signatures of the stau NLSP scenario, 
one also gets £ + £~ pairs or £ R tracks or decay kinks. 

Neutralino-stau co-NLSP scenario: If \rrif 1 — m^l < m T and < Tng R , both 

signatures of the neutralino NLSP and stau NLSP scenario are present at the same 
time, since Ni <->• fj decays are not allowed by phase space. 

Note that one always has rug > in the GMSB parameter space we explored, 
and that the classification we give above is only valid in the limit m e , m M — > and 
has to be intended as an indicative scheme. Indeed, we did not take into account very 
particular regions of the parameter space where, due to a fine-tuned choice of VF and 
the sparticle masses, one may achieve competition between phase-space suppressed decay 
channels from one ordinary sparticle to another and sparticle decays to the gravitino . 



Note also that we did not find in our sample any model with a sneutrino NLSP, since this 
is only possible in a corner of the parameter space where the lightest sparticle masses are 
well below 100 GeV. 

In Fig. [I], we show where in the (M mess , iV mess ) plane the scenarios described above are 
of relevance, for a GMSB spectrum of interest for the LC. One can see that the neutralino- 
NLSP scenario can only occur (but does not need to) for N mess = 1, 2 or 3. For iV mess = 3, 
it is also necessary to have a messenger scale as high as 10 8 GeV or more. Neutralino-stau 
co-NLSP models exist also for N mess = 4, but only for very high M mess . Stau NLSP and 
slepton co-NLSP models are instead possible for all allowed values of iV mess , but slepton co- 
NLSP models need M mess to be lower than 10 6 (10 7 ) GeV, if A mess = 1 (2). Perturbativity 
requirements up to the GUT scale start to be effective in excluding relatively low values 
of M mess for N mess > 5, while models with N mess = 7 or 8 are not possible if one imposes 
the simple cosmology-inspired condition itlq < 1 keV. 

Within a given scenario, the specific topology of the signatures is determined by the 
value of VF. We discuss this in detail in the next section, for the specific case where 
a neutralino is the NLSP. We now analyse a few important characteristics of neutralino 
NLSP models with 100 GeV < rrift < 250 GeV in our sample. 

First, in Fig. |^, we show the allowed regions in the (tan/?, A) plane for such a scenario 
to be realized. All the neutralino NLSP models we generated fall within the regions shown 
in Fig. ^| for a given value of A mess , but note that it is generally possible to construct models 
giving rise to different NLSP scenarios that also fall in the same regions of this plane. 
The regions in figure are sketched with a regular form to give a more intuitive feeling and 
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Figure 1: Scatter plot of GMSB models of interest for a 500 GeV e + e~ linear collider 
(100 GeV < mNLSP < 250 GeV) in the (M mess , N mess ) plane. The region in the upper left 
corner is excluded by requiring perturbativity up to the GUT scale. For each integer value 
of N mess (fractions along the y axis have here no meaning and are for display purposes 
only), from bottom to top and in different grey scale, we display neutralino-NLSP models 
(+), stau-NLSP models (*), slepton co-NLSP models (o), and neutralino-stau-coNLSP 
models (x). 
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Figure 2: Allowed regions in the (tan/3, A) plane for neutralino-NLSP models of interest 
for the linear collider, for different values ofN mess (necessary, but not sufficient conditions). 



are a bit wider than those actually populated by the relevant models in our sample. Also 
note that, due to the stau L R mixing producing lower mass eigenvalues for large tan/3, 
it is impossible to build a neutralino NLSP model of interest here for tan/3 > 30 (15), 
when the messenger sector is not the simplest possible one, namely N mess = 2 (3). 

Second, it is important to determine how much heavier the other sparticles can be 
compared to the N±. This tells us what the likelihood is that once A^-pairs are produced as 
an isolated signal at the LC, one can turn other SUSY processes on by just slightly raising 
the available c.o.m. energy. In all neutralino-NLSP models, the next-to-NLSP particles 
are the i?-sleptons, and in particular the fi (which turns always out to be dominated - 
85% or more - by the R component). The e R mass is particularly relevant, since it largely 
determines the e + e~ — > N\Nx cross section, together with the N% physical composition, 
due to the large contribution from t-channel e H -exchange graphs. Indeed, we will see that 
in GMSB models, the N% is dominated by the bino component, more strongly coupled 
to i?-particles, and the e L is always much heavier than the e R . The contribution from 
t-channel e L -exchange graphs is hence generally negligible. In Fig. |3], we show the ratio 
m SR /m^ i as a function of tan/3 and for different messenger multiplicity. We chose tan/3 
here as the independent variable mainly for visual purposes. The main information that 
can be extracted from Fig. |3| is that there are no neutralino-NLSP models where the e R 
is more than 1.8 (1.4, 1.2) times heavier than the N\ for A^ mess = 1 (2,3). 

Both in connection with the production cross section and with the decay properties to 
be discussed in the next section, it is essential to specify the possible physical composition 
of the Nx in GMSB models with neutralino NLSP. Fig. [| shows clearly that the bino 
component is always well above 90%, while the wino component never reaches the 2% 
level. The total higgsino component in the Ni can only rarely reach the 5% level, hence 
in some cases we will just neglect it in the following and assume that the Ni is a pure 
gaugino. Note also that in the EW-diagonalised basis, the photino component | (iVx | ^y) | 2 
of the NLSP is always included in the 0.60-0.85 range, while the zino component is in the 
0.15-0.35 range. 

The reader must be warned, however, that this is only true in the simple GMSB frame- 
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Figure 3: Bounds on the ratio between R-selectron and the N x masses in neutralino NLSP 
models of interest for the LC. Contours of the populated regions are shown for different 
values of iV mess . 
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Figure 4: Ni composition in GMSB models with neutralino NLSP of interest for the LC. 
From top to bottom and left to right, we show the bino, wino, higgsinoi and higgsino2 
components. Our basis for the neutralino mass matrix is the same as in Ref. 
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work we use in this paper. Here, we chose the phenomenological approach where one just 
assumes the existence of the \x and B terms at the messenger scale and determines them 
through EWSB conditions. As a consequence of this and the particular characteristics of 
the stop spectrum in GMSB, turns always out to be > 300 GeV 3> M z in neutralino 
GMSB models with 100 GeV <m^< 250 GeV. Further, the relation > 2M X always 
holds at the EW scale (here Mi is the bino soft mass). Such a circumstance produces 
decoupling between the gaugino and higgsino blocks in the neutralino and chargino mass 
matrices and the characteristic relations mg ~ m^ 2 ~ M 2 — 2m ^ ~ 2M\ approximately 
hold, while the heavier neutralino and chargino mass eigenvalues are always of order 

However, there are many possible sources of more complex scenarios. For instance, 
if one attempts to put together a radiative mechanism to generate and B, one may 
find extra corrections to the Higgs soft (mass) 2 parameters, which in turn can modify the 
value of \fi\j often lowering it |8], |9|, [□], p2[ . Further, it is possible to build coherent GMSB 
models that have unequal messenger multiplicity relative to the three SM gauge groups. 
Models in this class exist where the higgsino component of the Ni NLSP is large, with 
remarkable phenomenological consequences, both for the NiNi production cross section 
and Ni decay BR's to be discussed below. Many other variations in the messenger sector 
are possible |23| , but the associated phenomenology is beyond the scope of this paper. In 
the following, we will always assume that the higgsino components of the N\ NLSP are 
small and possibly negligible. 

We list here three reference GMSB models with neutralino NLSP in our sample that we 
will use in the following. We chose these particular models because they are qualitatively 
different for our experimental studies of Sec. [7|. They cover a good spectrum of possibilities 
and provide a feeling of the various problems that the experimenters could face if nature 
had chosen GMSB and the neutralino as the NLSP. 

Model # 1 features are summarised in Tab. [3], where spectrum, production cross 
sections at a 500 GeV LC, and other relevant details such as sparticle physical composition 
and dominant decay channels are given (details about the NLSP decay are deferred to the 
next section). The precise values reported for masses and cross sections (which include 
ISR and running a em effects), depend slightly on details of the spectrum calculation, 
higher-order corrections etc. They should be considered as an approximation at the level 
of a few percent. (Since here we are not particularly interested in the Higgs sector, we 
give only indicative information about it.) 

Model # 1 is a model with a rather light spectrum, in particular the NLSP mass is right 
at our assumed LEP2/Tevatron bound of 100 GeV. The next-to-NLSP, the i?-sleptons, are 
in the middle of their allowed mass range for such a NLSP mass. If this GMSB scenario 
were to be realized, the sparticles that could be produced with appreciable cross section 
at a 500 GeV LC would be gauginos, i?-sleptons and L-selectron (in association with 
e R ). The total GMSB signal would be in this case quite "generous." Heavy interacting 
sparticles are definitely out of reach, even for a possible second phase of LC operations 
with c.o.m. energy at or slightly above 1 TeV. These large mass splittings are a well- 



known characteristics of GMSB models (cfr. e.g. Refs. || |24j) and are due to the fact 
that gaugino and scalar masses are proportional to the relevant gauge couplings. (The 
light Higgs is close to the edge of detectability at LEP2/Tevatron, depending on fine 
details and higher-order corrections that we do not take into account. In any case, models 
with a slightly heavier h° and no significant differences in the other sectors can easily be 
constructed with small changes to the input parameters. The rest of the Higgs sector is 
very heavy.) Notice that such a model would not be expected to produce a large signal 
at the LHC, due to the heaviness of gluino and squarks, hence a careful search and study 
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Model # 1. INPUT: M mess = 161 TeV; N mess = 1; A = 76 TeV; tan/3 = 3.5; // > 


Particle 


Mass 


Production @ 500 GeV LC 


Comments 


G= LSP 


> 2.9 eV 


indirect only 


stable 


Ni = NLSP 


100.0 GeV 


<7(JViiVi) = 256 fb 


|(JVi|B) 2 = 0.97; decays to G 




136.6 GeV 


a(f l f l ) = 56.9 fb 


~ f R ; decays to Nx 


e R , fiR 


137.1 GeV 


a(e R e R , fl R fi R ) = 274,56.6 fb 


decay to Nx 


C x 


183.3 GeV 


aidd) = 137 fb 


\Uxi 


| 2 = 0.87; |y n | 2 = 0.94 


N 2 


184.6 GeV 


a(NxN 2 , N 2 N 2 ) = 39.1,38.3 fb 


|(Ar 2 |W 3 )| 2 = 0.9 




264.4 GeV 


— — 


v T slightly lighter 




274.3 GeV 


a(e L e R ) = 101 fb 




^2 


274.5 GeV 


a{fxf 2 ) < 0.1 fb 


~ f L 




~ 105 GeV 


a(hZ) ~ 70 fb 




iV 3 , A> 4 , C 2 


> 400 GeV 






#°, v4°, #± 


> 500 GeV 






9 


~ 650 GeV 






Q 


> 700 GeV 







Table 1: Input parameters, output spectrum and basic characteristics of a typical GMSB 
model with a 100 GeV neutralino NLSP: Model # 1. 

at the LC would be most likely necessary, if not for initial SUSY discovery, then at least 
for a confirmation and for determining with good accuracy the source of the anomalous 
signal and the underlying SUSY-model parameters. 

Model # 2 (see Tab. 0) is much more of an "avaricious" model, with a 200 GeV NLSP 
mass. It is obtained from Model # 1 by just raising the input value of A, leaving the 
^mess/A ratio and the other parameters untouched. The only GMSB signal present at a 
500 GeV LC would be NLSP pair production in this case. Note that changing the iVi mass 
with respect to Model # 1 does not only result in a drastic reduction of the cross section 
for NiNi production, but also in an important change in the N± decay BR's, as described 
in detailed in Sec. || As a consequence, even focussing on A r iA r i production only, this 
model would produce a considerably different signal at the LC compared to Model # 1, 
both quantitatively and qualitatively. In this case, gluino and squarks are very heavy, 
possibly close to a reasonable bound from naturalness arguments and the GMSB signal 
at the LHC would be rather scarce. 

Model # 3 (see Tab. is a special model presenting some unusual and challenging 
characteristics. First of all, the i?-slepton masses are very close to the neutralino NLSP 
mass of 165 GeV. (Note that when the difference between the f\ mass and the NLSP mass 
approaches the tau mass, one falls in the neutralino-stau co-NLSP scenario that we are 
not treating here.) As we will see, this poses the problem of separating the various GMSB 
signals from each other at the LC in order to perform specific measurements. Another 
experimental challenge follows from the fact that the relatively low minimum gravitino 
mass combined with a quite large N\ mass makes it possible for the neutralino to decay 
very close to the interaction region at the LC (see Sec. ^) in this case. 

In Fig. [|, we plot the cross-section for the various SUSY production processes as a 
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Model # 2. INPUT: M mess = 309 TeV; N mcss = 1; A = 146 TeV; tan/3 = 3.5; fi > 


Particle 


Mass 


Production @ 500 GeV LC 


Comments 


G= LSP 


>H eV 


indirect only 


stable 


iVi = NLSP 


200.0 GeV 


<t(JViJVi) = 42.3 fb 


{(N^B) 2 = 0.99; decays to G 


h° 


~ 115 GeV 


a(hZ) ~ 63 fb 




fi 


256.4 GeV 


— 


~ tr\ decays to N\ 




256.8 GeV 


— 


decay to Ni 


c, 


374.1 GeV 


— 


C/nl 2 = 0.95; \/ n 2 = 0.98 


N 2 


374.4 GeV 


— 


|(jV 2 |iy 3 )| 2 = 0.96 




511.5 GeV 


— 


v T slightly lighter 


e L , P>l 


516.7 GeV 


— 






516.7 GeV 




- r L 


N 3 , N 4 , C 2 


> 700 GeV 






H°, A , H± 


> 900 GeV 






9 


~ 1150 GeV 






Q 


> 1300 GeV 







Table 2: Input parameters, output spectrum and basic characteristics of a typical GMSB 
model with a 200 GeV neutralino NLSP: Model # 2. 



Model # 3. INPUT: M mess = 110 TeV; iV mcss = 1; A = 100 TeV; tan/5 = 3; // < 


Particle 


Mass 


Production @ 500 GeV LC 


Comments 


G= LSP 


> 2.6 eV 


indirect only 


stable 


^ = NLSP 


165.0 GeV 


<r(JViiVi) = 136 fb 


(A/"iLB) 2 = 0.99; decays to G 




171.5 GeV 


<7(fifi) = 34.6 fb 


~ f^; decays to Ni 


e R , J1r 


171.8 GeV 


<r(e R e R , PlrPlr) = 78.7,34.5 fb 


decay to Ni 


N 2 


315.0 GeV 


(j(A>iA>2) = 1.07 fb 


\(N 2 \W 3 ) 2 = 0.96 


h° 


~ 105 GeV 


a(hZ) ~ 70 fb 




Ci 


315.1 GeV 




f/ n 2 = 0.93; |y n | 2 = 0.99 




342.5 GeV 




v T slightly lighter 


e L , Jil 


349.8 GeV 








349.8 GeV 






N 3 , A> 4 , C 2 


> 500 GeV 






H°, A , H± 


> 650 GeV 






9 


~ 950 GeV 






Q 


> 950 GeV 







Table 3: Input parameters, output spectrum and basic characteristics of a typical GMSB 
model with a very-short lived neutralino NLSP and nearly degenerate light sparticles: 
Model # 3. 
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Figure 5: Inclusive signal "7" "7" + X + $ from GMSB Model # 1 as a function of y/s 
in the range of interest for a 500 GeV LC. The total as well as all contributions to the 
signal from each sparticle-pair production process are shown. The normalisation is based 
on an integrated luminosity of 100 fb _1 . 

function of y/s in the range of interest for a 500 GeV LC, for Model # 1. Actually, Fig. [5] 
contains some more information, since the normalisation of the y axis takes into account 
the inclusive nature of the GMSB signal and the typical luminosities of the LC. As we 
will see in Sec. |3], for all neutralino NLSP models N\ — > G7 is the dominant NLSP decay 
channel, with BR's always greater than about 85%. (For Model # 1 this is actually 95%.) 
As a consequence, each time a sparticle pair is produced, there is a large probability of 
getting a final state with 2 photons, some other particle resulting from cascade decays 
and large missing energy. In Fig. [5], the quotation marks for 7 mean that (one of) the 
photons might escape detection if the N\ lifetime is very large (cfr. Sec. |3]). The number 
of "7" "7" + X + Jfi events is normalised to an integrated luminosity of 100 fb _1 , which 
is a typical value for a few-month run at the LC (cfr. Sec. |4j). Notice here that for the 
case of Model # 1, running at c.o.m. energies of order 270 GeV would still allow for order 
20,000 GMSB events, while selecting pure 77-^ events only. This circumstance will be 
exploited for our experimental studies in Sec. 0. We will also use Fig. |] as a basis for the 
study of Sec. [|. 

3 Neutralino NLSP Decays 

In this section, we analyse the properties of the NLSP decay in GMSB models, with focus 
on the case where Ni = NLSP. 

In Ref. [|IIJ, all the formulas for 2-body decays involving the gravitino can be found, in 
the limit where the gravitino interactions can be approximated by those of the goldstino 
and its mass can be kinematically neglected, which is always the case in GMSB at collider 
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Decay Channel 
N t -> SG 


A s = 


where 


S = 1 




K i7 = \Nn cos 9 W + N i2 sin 9 W 


2 


z° 




^iz T = \Nn sin 9 W - N i2 cos 9 n 
Kiz L = \N i3 cos/3 - Na sin (3\ 


|2 

2 


h° 




Kih° — 


N i3 sin a — Nu cos a 


2 


H° 


K iH o/2 


KiHO = 


N i3 cos o; + Na sin a 


2 


A 


n iA o/2 


KiAO = 


A" i3 sin/5 + iYj 4 cos/? | 2 



Table 4: Constant factors entering the expressions for the widths of neutralino 2-body 
decays to gravitino. 



energies. For a generic decay S — > SG, where S is a SM particle and S its MSSM 
superpartner, one has for the corresponding width, 

r -— 4i(® 8 44(^ (4) 



48tt M' P 2 ml ■ b ' 16ti y/p 



where the gravitino mass is given by (]I|), f3* s is the relativistic factor Jl — (ms/m§) 2 if S 
is a vector or scalar boson. If S is a massless fermion, — ► 1. ^4^ is a constant depending 
on the S, S spin and possibly a mixing matrix element. For example, if S is a SM lepton 
or quark and S a slepton or squark, then simply As = 1. 

We are here interested in the S = N\ case, since in our neutralino NLSP models the 
only particle that can undergo a 2-body decay to a gravitino with a non-negligible width 
is the lightest neutralino. The relevant expressions for As can be found in Tab. |1 where 



we use the notation of Ref. |21J for the neutralino mixing matrix and a is the mixing 
angle in the MSSM neutral Higgs sector. 

Due to the absence of the (/3*) 8 kinematic suppression and the N\ physical composition 
in the models of interest here (cfr. Fig. f|), the Ni decay is always dominated by the photon 
channel. However, in the context of this paper where > 100 GeV and fine details 
of the neutralino decay will be used in the following, it is important to note that the BR 
for decaying to a Z° can be sizeable and also to address the problem of the 3-body decay 
channels Ni — > f fG, where / is a SM lepton or quark 0. 

The Feynman diagrams contributing to the 3-body processes are shown in Fig. || An 
analytical formula for the sum over final-state SM-fermions of the total widths for these 
decays via real or virtual boson exchange has appeared in Ref. . However, that summed 
formula could not take into account Diags. 6-9 where a L- or i?-sfermion is exchanged, 
which are a-priori not less relevant than Diags. 2-5, where other heavy intermediate 
particles are involved. Also, in Ref. the virtual photon contribution (Diagr. 1) was 
calculated including an overall detector-dependent cutoff on the fermion pair invariant 
mass, chosen to be of order 1 GeV. In the context of this paper, however, we have at our 
disposal a full detector simulator (cfr. Sec. ^) and we will be interested in the individual 
BR's for each // pair (cfr. Sec. 0). Also, the kinematical distributions of these decays 
are relevant to our following studies. In the rest of the paper, we will often use the name 
"neutralino charged decays" when referring to the N\ — > f fG channels and in particular 
to the case where the subsequent final state includes either a charged lepton or a charged 
"stable" hadron. 
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Diagr. 1 



Diagr. 2 



Diagr, 3 



Diagr. 4 



Diagr. 5 



Diagr. 6 



a 



N , 



fit i^f 



Diagr. 7 



Diagr. 9 



Figure 6: Feynman diagrams contributing to 3-body neutralino decays Ni — > ff'G, where 
f is a charged SM fermion. The channels Ni — > VjVjG receive contribution from Diags. 2, 
6 and 7 only. 



To account for these channels, we proceeded as follows. Using the general lagrangian 
for the goldstino interactions with matter without assuming on-shell conditions (cfr., 
10 1), we input all the relevant vertices^ involving the gravitino in CompHEP 



e.g. 
3.3 



Refs. 
18 [Ei 



in a limit suitable for collider physics. For the other vertices involving MSSM 
particles, we used the home-made lagrangian^ that was first checked against analytical 
calculations and then used for numerical evaluations in the work of Ref. [ 20 1 

We named the resulting software Gravi-CompHEP [|l]|]. Using Gravi-CompHEP, we 
found that the contribution to the total width from virtual photons is in very good 
numerical agreement (for the electron and muon case up to 4 digits) with the analytical 
formula [28], valid to lowest order in mf/m^ 



r(Ni - ffG) = r(iv 1 



jG) — 
Sir 



N fQ} 



In 



m N 1 
. m f 



15 
T 



(5) 



where Qf is the final fermion electric charge in units of e, Nj = 1(3) for leptons (quarks) 
and the cutoff is naturally provided by the / mass. Note that, e.g. for a 100 GeV 
neutralino mass as in Model # 1, Eq. (|J) gives for the case of electrons in the final 
state numbers about twice (five times) as large as for the case of muons (taus). For 
hadronic final states, a realistic evaluation must take hadronization effects and higher- 
order corrections into account. However, since hadrons will not be our main focus in the 
analyses of Sec. |7| and we will be most interested in the BR's for the leptonic channels, 



^or some of our gravitino vertices and using the goldstino lagrangian as an input, we checked that 
there is agreement with the output of LanHEP 1 . 5 . 06 ]25| ] 

2 For the relevant MSSM vertices and in the relevant limit, we checked that there is numerical agreement 
with our results when using the CompHEP 3 . 3-compatible MSSM lagrangian of Ref. (27) instead. 
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we chose a reasonable approximation using a rough cutoff for the invariant mass of the 
final fermion pair at 2Aqcd ~ 300 MeV for light quarks and Eq. (H) for heavy quarks. 

The contribution from Z°-exchange to Y(N\ —>■ ff'G) (Diagr. 2) is obtained from 
Eq. (f|) by replacing 



BR(Z-//), 



where Jo and l\ are kinematical factors taking finite Z°-width effects into account^]. In 
our model sample for the LC, one finds that as long as > 120 GeV the on-shell Z° 
approximation is accurate at the level of 10% or better. For lighter neutralinos, the full 
calculation is required, since e.g. for = 100 GeV one has Iq(Z) = 0.0052 and h(Z) = 
0.0023, whereas (/3^) 8 = 0.00081, and the on-shell Z° approximation underestimates the 
Diagr. 2 contribution by a factor 2.5-3. On the other hand, for models with rrifj < 120 
GeV in our sample, the Z°-exchange contribution is always < 3 (15)% of the virtual 
photon contribution, e.g. for the / = e (r) case. 

The interference between the 7- and Z°-exchange diagrams turns out to be always 
small, generally at the level of a few % or less of the pure Diagr. 1 contribution (cfr. also 
Ret 1). 

Diagr. 3 is always negligible for the models of our interest here. Indeed, one has both a 
dynamical suppression due to the lack of higgsino components in the N\ (cfr. Fig. f|) and 
a kinematical ~ (/5*) 8 suppression, for ^ 150 GeV (in this range, nth/m^ < 0.85 
always in our model sample and the on-shell approximation applies). Taking off-shell 
effects into account for < 150 GeV (the formulas are similar to those for the Z° case 
described above) does not help either, since the h° width is typically very small in the 
MSSM. Also, Diagr. 3 contributes to the channels with heavy fermions in the final state 
only, which have typically lower BR's. Diags. 4 and 5 are even more strongly suppressed, 
because the masses of the CP-odd and heavy CP-even Higgses only rarely drop below 
300 GeV in our model sample. Interferences involving Diags. 3-5 are basically zero. In 
the rest of this section, we will often assume for simplicity that the N\ is pure bino; this 
makes all the contributions from Diags. 3-5 zero and is justified by Fig. [|. 

Finally, as far as Diags. 6-9 are concerned, the / exchanged is necessarily heavier 
than the initial N\ by definition of neutralino NLSP model. For the case of hadronic 
final states, these diagrams do not count, since the squarks are too heavy in our models. 
However, it turns out that limited to the case of Ir exchange, the contribution to the 
width is often non-negligible and at the level of several to 10% of the total, especially for 
those models where m^Jm~ lR is close to 1 and for the case of heavier leptons in the final 
state where Diagr. 1 is less dominant. This is again due to the relatively large B — £r 
coupling and the fact that the £r can never be much heavier than the NLSP (cfr. Fig. |3[). 
Diags. 6-7 are always negligible, even in the leptonic case, due to the relative heaviness 
of the L-sleptons. 

In Fig. |7|, we give a general idea of the behaviour of the BR's for the main neutralino 
NLSP decay channels as a function of for all the models in our sample of interest for 
the LC. From top to bottom, we show the BR's for the dominant two-body channel to a 
photon, the two-body channel to a Z° including off-shell effects (so that the contribution 
from Diagr. 2 to 3-body channels can be readily extracted by multiplying by the appro- 
priate Z° BR), the hadronic and e + e~ 3-body channels from Diagr. 1. For comparison, 
we also report our results for the BR of the 2-body Ni — > h°G decay in the on-shell 



3 Analytical expressions can be found in Refs. |], 
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approximation. The logarithmic scale does not allow inspection of fine effects. However, 
it is evident that the Z° channel can be important with BR's up to about 15% for heavy 
neutralinos, while the main 3-body channels via virtual photon are always relevant with 
BR's at the level of a few %. The Higgs channel has always BR's less than 0.1%. Note 
also that, due to the homogeneous physical composition of the Ni in our sample, the 
important BR's are basically dependent only on the neutralino mass. 

Fig. P is a scatter plot for our model sample showing in detail the BR's for all the 




Figure 7: Scatter plot for the BR 's of various neutralino NLSP decay channels as a function 
of the Ni mass. Dots in different grey scale (colours) refer to the decays Ni — > 7G, 
N\ — ► ZG (including off-shell effects), and to hadrons or e + e" via virtual photon, as 
labelled. For reference, we also report results for the two body Ni — > h°G decay in the 
on-shell approximation, whose BR is always negligible. 

3-body Ni decay channels (excluding Ni — > vvG) as a function of the neutralino mass. 
From top to bottom, hadronic + G, e + e~~G, fi + fi~G and t + t~G final states are calculated 
including all contributions from Diags. 1-9 in Fig. The BR's for all 3-body channels 
all increase for heavier neutralinos. The hadronic channel occurs about 2% to 15% of 
the times, the electron channel 1%— 1.5%, the muon channel 0.8%-l.l%, the tau channel 
0.3%-0.7%. Again, fixing the N\ mass basically determines these BR's, with some more 
uncertainty for the hadron and r channels that receive relatively larger contributions from 
Z°-exchange. 

In spite of the fact that the BR for the leptonic 3-body channels is often quite low, due 
to the large integrated luminosity that might be available at a LC, the possible number 
of events featuring a (displaced) £ + £~ pair is still large in most cases. In Fig. |9|, we show 
scatter plots for our neutralino NLSP model sample referring to inclusive GMSB signals 
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Figure 8: Scatter plot for the BR's of neutralino NLSP decays to several 3-body ffG 
channels as a function of the Ni mass. Dots in different grey scale (colours) refer to 
decays to quarks, electrons muons and taus, from top to bottom. Contributions from all 
Diags. 1-9 in Fig. |6| are included here. 



at a y/s = 500 GeV LC as functions of the neutralino mass. We refer to a nominal 100 
fb _1 run and sum over all SUSY production processes. In Fig. |9](a) we report the number 
of events including two (displaced) photons and missing energy coming from two long- 
lived neutralino decays. In Fig. |9|(b), we consider all events including at least a displaced 
e + e~ pair and missing energy. In the models we are interested in, one gets at least 100 
such events if rrift < 200 GeV, and up to about 10,000 events for lighter neutralinos. 
Notice that the meaning of "displaced" here is that the tagged particles are produced at 
some distance from the interaction region, where the distance depends on the neutralino 
lifetime and the specific processes considered. In some cases, the displacement might be 
so large that the particles are actually produced outside the detector, as we will see in 
the following. 

In Sec. 0, we will see that it is sometimes useful to run a LC at a c.o.m. energy not 
far from the NiNi threshold to isolate the signal from neutralino pair production. To 
give a feeling about this problem, in Fig. [lOj we show a scatter plot similar to Fig. 
but for y/s = y/2rrif ]i . The big black dots refer to models for which only N\Ni pairs can 
indeed be produced at such an energy, while small grey dots are for models where other 
processes (typically pair production of i?-sleptons) are also below threshold. Note that 
the number of events including a displaced e + e" pair is always larger than about 100 for 
an integrated luminosity of 100 fb -1 . The BR's for all the main N± decay channels for 
the three reference models introduced in Sec. |2| are shown in Tab. [| and will be referred 
to in the analyses of Sec. 0. 
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Figure 9: Scatter plot for some inclusive GMSB signals in neutralino NLSP models of 
interest for a LC. We report the number of events expected after a nominal 100 fb^ 1 
run at 500 GeV c.o.m. energy, (a) Events including two displaced photons and missing 
energy, (b) Events including a displaced e + e~ pair and missing energy. Here we assume 
that the delayed decays of the N% NLSP all occur within the detector. 



Decay Channel 


BR in Model # 1 


BR in Model # 2 


BR in Model # 3 


Ni 


0.9507 


0.8395 


0.8913 


Nx -> ZG (a) 


0.0003 


0.1115 


0.0585 


Ni -> e + e~G ^ 


0.0164 


0.0191 


0.0179 


Ni -> n+n~G < b > 


0.0079 


0.0117 


0.0101 


iVi -> r+r-G ( b ) 


0.0034 


0.0076 


0.0059 




0.0213 


0.0999 


0.0631 


Nx^YtWiG {b) 


0.0002 


0.0223 


0.0117 


iVi -> h°G ^ 




< 0.0001 


< 0.0001 



Table 5: BR's for the main decay N\ decay channels in our three reference GMSB models. 
( a ) Entries include the on-shell Z° (h°) contribution only. ( b ) Entries include all contribu- 
tions. As a consequence, due to double counting of the on-shell Z° (h°) contributions, the 
BR's in each column do not sum up to 1. 
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Figure 10: Scatter plot for some inclusive GMSB signals in neutralino NLSP models of 
interest for a LC, as in Fig. |], but for yfs = y^m^. Big black dots are for models where 
only N]_Ni pairs can be produced at such an energy, while little grey dots refer to models 
where other processes are also below threshold. 



20 S. Ambrosanio, G. A. Blair / Measuring GMSB Parameters... 




M(i + n/M N1 M(i + n/M, 



Figure 11: Normalised £ + £~ invariant mass distributions for the leptonic three-body de- 
cays of the Ni in Model # 1 (a) and Model # 2 (b). Stars, circles, crosses refer to the 
electron, muon, and tau case, respectively. 

In Sec. 0, we will heavily use the characteristics of the three-body decays of the neu- 
tralino, including their kinematical distributions. Using Gravi-CompHEP, we calculated 
such distributions for our reference models and we show here our results. We then im- 
plemented numerically the corresponding differential widths into our event generator to 
perform the Monte Carlo simulation. 

In Fig. [11], the normalised invariant £ + £~ mass distribution for all the leptonic channels 
is plotted. The stars, circles, crosses are the central points of our results for the electron, 
muon, tau cases respectively in Model # 1 (a) and Model # 2 (b). The horizontal bars 
show the Mn/ m m binning we used, while the errors coming from numerical phase space 
integration on the distributions are too small to be visible in logarithmic scale in most 
cases. Note that the distributions are sharply peaked for low invariant masses close to 
2Mu and this is more and more true for lighter leptons. As a consequence, e.g. the e + e~ 
pairs coming from NiNi production at the LC and subsequent three-body decay of (one 
of) the neutralinos tend to be generated with small separation angles, which introduces 
some experimental challenges (cfr. Sec. [F]). In Fig. |TT|(b), the peak corresponding to 
the Z°-exchange contribution is also evident. To allow a better inspection of the scaling 
properties of the distributions with and mj, we used here the (very good) N\ = B 
approximation for both models, so to keep the N\ physical composition constant when 
going from Model $ 1 to Model # 2. For the case of Model # 1, we included Diag. 1 
only of Fig. || since the other contributions would hardly be visible in the plot anyway. 

Limited to the case of Model # 1, in Fig. O, we show some relevant angular distri- 
butions for the Ni — > e + e~G decay, the three-body channel we will be most interested in. 
The circles refer to the normalised cos8(e + e~) distribution, where 8(e + e~) is the angle be- 
tween the electron and the positron momenta in the decaying Ni rest frame. As expected, 
the e + and the e~ prefer to proceed along the same direction, due to the dominance of 



S. Ambrosanio, G. A. Blair / Measuring GMSB Parameters. 



21 



a 




* , - 



-0.25 0.25 0.5 0.75 

cos8,6 (e"V,eG) 



Figure 12: Normalised angular distributions for the N± — > e + e decay in Model # 1. 
Circles, squares, crosses and stars refer to different angles, as defined in the text. 



the virtual photon contribution. The stars correspond to the normalised cos8(e ± G) dis- 
tribution, where 8(e ± G) is the angle between the electron (or positron) and the gravitino 
momenta in the Ni rest frame. The e + e~ pair is produced in the direction opposite to 
the G in the great majority of cases. The squares refer to the normalised cos#*(e + e") 
distribution, where 8*(e + e~) is the angle between the electron (or positron) momentum 
and the direction of the boost of the e + e" system with respect to the Ni rest frame, 
calculated in the e + e~ rest frame. In our case, this is basically the angle between the elec- 
tron (or positron) and the virtual photon momenta and the almost constant behaviour is 
then expected. Finally, the crosses show the cos8*(e ± G) normalised distribution, where 
8*(e ± G) is defined as above and refers to the e ± G system instead of the e + e" one. The 
general behaviour of these angular distributions is better understood in the light of the 
fact that the two-body Ni — > 7G decay is isotropic. Some instability of our results due to 
numerical phase space integration is visible, but it is well within the shown vertical error 
bars. 

In Fig. [13], we show the normalised energy distributions for the Ni — > e + e~G decay in 
Model # 1 (a) and Model # 2 (b). The stars refer to the electron (or positron) energy, 
while the circles are for the gravitino energy, in the decaying N\ rest frame. For the 
case of Model # 1, the G tends to take half of the available energy, while the electron 
and positron share the rest with an almost uniform distribution between m e and about 
Model 7^ 2 features evident effects of the Z°-exchange contribution that add to a 
behaviour similar to the one for Model # 1 coming from the still dominant virtual photon 
contribution. Again, in order to allow a cleaner comparison between the two models, 
we used here the bino approximation. Note that both the scales on the x and y axes in 
Fig. [13] are interrupted for display convenience. 

After having inspected the various possible decay channels, we turn now to the dis- 
cussion of the total width and the lifetime of the neutralino. As anticipated above, this 
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Figure 13: Normalised electron or positron (stars) and gravitino (circles) energy distribu- 
tions for the N~i — > e + e" decay in Model # 1 (a) and Model # 2 (b). 



determines the topology of the signatures of neutralino GMSB models at colliders. A 
single neutralino produced with energy will decay before travelling a distance A with 
a probability given by 

F(A) = l-exp(-A/L) where (6) 
L = cr^ipj)^. (7) 

L is the Ni "average" decay length and (f3y) ^ is the kinematical factor (E 2 ^ /m 2 ^ — l) 1 / 2 . 

Note that for Ni pairs directly produced at the LC with y/s = 500 GeV, {Pyjfa = 2.29 
(0.75) if = 100 (200) GeV. 

Using Eqs. (|l|), (f|), and (||), the neutralino lifetime can be conveniently expressed in 
the suggestive form, 

16vr VF* 

1 ( \*( m r 5 (8) 

100i3 V100 TeV; V 100 GeV/ ' K} 

which stresses the scaling properties with the 5 th inverse power of the neutralino mass and 
the 4 th power of the fundamental SSB scale. B is a number of order unity that can be well 
approximated by B ~ = when the two-body Ni — ► 7G channel widely dominates, 
as e.g. in Model # 1, or by simple expressions in most cases. In general, however, it is 
a complicated function of the neutralino composition, the GMSB model spectrum etc., 
when the full contributions to the three-body channels are taken into account. 

Once the neutralino mass and lifetime are measured (cfr. Sec. |^), one can get striking 
information on \fF from Eq. (§). The uncertainty is then only due to the factor B. If 
the BR's for the various N\ decay channels are also measured with good precision, or the 



CT 



S. Ambrosanio, G. A. Blair / Measuring GMSB Parameters... 



23 



torfrafc© UL8P USMIam -» V# 



(a) 



#4 




• •• m(e R )<150GeV 
m(e R )>150GeV 



UL O.I 



120 140 160 180 200 220 240 

mP,) [GeV] 



(b) 




• •• m(N 1 )<150GeV 



mfl^) > 150 GeV 



200 250 300 350 400 

m(i R ) [GeV] 



Figure 14: Scatter plot showing the relation between the neutralino NLSP lifetime and 
the fundamental scale of SUSY breaking \/~F , i.e. the factor B in Eq.(jj$), as a function of 
the neutralino mass (a) or the R-selectron mass (b) in GMSB models for the LC. Big grey 
dots in (a) represent neutralino NLSP models with a light R-selectron (102-150 GeV), 
small black dots are for the heavier selectron case (150-430 GeV). In (b), big grey (small 
black) dots are for models with mfj < (>) 150 GeV. 



neutralino composition and the (light) GMSB spectrum is extracted by measuring other 
observables (production cross sections, distributions), the fundamental SUSY breaking 
scale can be determined precisely. However, it is remarkable that even without collecting 
additional information, the knowledge of and cr^ is sufficient to constrain the value 
of yF in a narrow range, based on the well defined characteristics of GMSB models. In 
Fig. 0, we report scatter plots of our neutralino NLSP model sample for the LC showing 
B in Eq. as a function of the neutralino mass (a) and the right selectron mass (b). To 
stress the existence of some correlation between B and the right selectron mass for models 
with a fixed light neutralino mass, in Fig. |l4](a), the big grey (small black) dots refer to 
models where 102 (150) < rrig R < 150 (430) GeV. This can be compared to Fig. |T4|(b), 
where big grey (small black) dots correspond to models with < (>)150 GeV. From 
this, one can e.g. infer that if — 120 GeV and the neutralino lifetime is measured to 
be about 1 cm, then 360 < yF < 385 TeV. If, in addition, m e - fl is measured to be heavier 
than about 150 GeV (for instance, from e R e R threshold scanning, see Sec. || or using its 
impact on the A^A^ cross section), then the allowed range is further reduced to 370-385 
TeV. For a 200 GeV neutralino, a 1 cm lifetime gives 725 < v"F < 740 TeV. 

To summarise, we note that in the absence of further information, the theoretical error 
on determining \fF from given values of N\ mass and lifetime amounts to about 3% in 
the worst case, helped by the 4 th power dependence in Eq. (§). 

It is of primary importance for collider phenomenology to assess the range of variation 
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Figure 15: Scatter plot of the neutralino NLSP lifetime as a function of the messenger 
scale M mess (a) and (b). For each set of GMSB model input parameters (A, M mess , 
etc.), the lower limit of the NLSP lifetime is plotted, corresponding to y/F ~ V^mess — 
\J AM moss . Only models that fulfil the limit on the gravitino mass (mg < 1 keV =>- 
V^mess ^ VF < 2000 TeV) suggested by simple cosmology are used. 



for cr^ . As anticipated in Sec. |], it is possible to use a lower limit from theory on y/~F, 
while significant upper limits can only come from weak cosmological arguments suggesting 
rug, < 1 keV. The lower limit defines a minimum value cr^ in for the neutralino lifetime 
as well as m^ m for the gravitino mass on a GMSB model- by-model basis. In Fig. 15, we 
plot this limit as a function of M mess (a) and (b) for our model sample of interest 
for the LC. In this plot, we also use the cosmological upper limit on mg, in the sense 
that those models where m% m > 1 keV are not plotted. As a result, one can see that the 
neutralino lifetime can be anywhere between about 5 microns and about 25 metres (or 
more if no cosmological arguments are used). Models with a high messenger scale produce 
longer neutralino lifetimes. For instance, if M mess > 10 4 TeV, then cr^ is always larger 
than about 1 cm. Also note that shorter lifetimes are obtained for heavier neutralinos. A 
100 GeV neutralino will always live more than about 15 microns. On the other hand, a 
250 GeV neutralino will tend to decay well within a typical detector size, with lifetimes 
always smaller than about 20 cm, if cosmological arguments are used. 

When SUSY pairs are produced in a neutralino NLSP scenario, the resulting final 
states always include two neutralinos, which in turn decay to a gravitino + X. The prob- 
abilities of both, one or zero neutralinos decaying within a given volume of the detector 
depend on the neutralino's decay length L of Eq. (0), which in turn depends on the spe- 
cific SUSY process, model and collider c.o.m. energy one is considering. If we define a 
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Figure 16: Number of events in 100 fb featuring one or both neutralinos decaying within 
two reference spheres of radius 1.6 and 3.0 m (different grey scale or colours) after NiNi 
production at y/s = 270 GeV as a function of L or y/F for Model # 1. 



spherical volume of radius R, then the probabilities associated with these circumstances 
are of course given by P(R) 2 , 2P(R)[1 - P(R)], and [1 - P(R)} 2 . In Fig. [L|, we 



show how many events are expected with two (solid line) or one (dashed line) neutralino 
decays as a function of L (or y/~F) for the case of direct neutralino-pair production in 
Model # 1. We show curves for two reference spheres with a radius of 160 and 300 cm 
(we will see in Sec. ^| that the outer cylinder of a typical proposed TPC for a LC detector 
is included between such two spheres). Our numbers refer to a 100 fb -1 run at the LC 
with y/s = 270 GeV, where only NiNi pairs can be produced with a cross section of 
188 fb. This choice of parameters will be of relevance for the studies to be presented in 
Sec. 0. Note that for neutralino decay lengths as large as 1 km (and \J~F larger than 5000 
TeV), there still are about 100 events where one neutralino decays within the reference 
volume. We warn however that this is only true if a sum over all possible final states 
coming from neutralino decays is performed and no angular or other detector acceptance 
cuts are taken into account. Based on this, a refined statistical study based on a realistic 
cylindrical detector and experimental framework for the LC is performed in Sec. [7.7| . 
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4 The Linear Collider and the TESLA Project 

The LHC will explore the next high energy frontier and can be expected to be among 
the prime sources of new physics discoveries into the next decade and beyond. However, 
ongoing studies of the physics potential of a e + e~ LC operating at c.o.m. energies ranging 
up to 500 GeV, or higher (1-2 TeV), are revealing many complementary measurements 
that could be made at such a machine on a similar timescale to that of the LHC. Addi- 
tional options available at a LC are a considerable electron (and possibly also positron) 
polarisation, e^ and 77 options, as well as the potential for e~e~ collisions, making a LC 
a very flexible and relevant facility, with particular application to detailed studies of new 
physics signals. 



Several linear collider designs are presently under discussion [2!J and much of the dis- 
cussion in this paper is applicable to any machine. However, in order to relate our study to 
a specific case, we explore the machine parameters of the high-luminosity TESLA option. 
The TESLA machine proposal is described in some detail in vol. 2 of the ECFA/DESY 
"Conceptual Design Report" (CDR) [P0[| . The most recent proposals involve two phases of 
operation; an earlier phase operating at y/s = 500 GeV or less and a later phase operating 
at y/s = 800 GeV or less, with luminosities of 3.1 x 10 34 cm _2 s _1 and 5.0 x 10 34 cm^ 2 s _1 
respectively [pl|. In this study, we mostly limit ourselves to the first phase foreseen for 
such a collider, with y^i varying between approximately 200 and 500 GeV. In most cases, 
we will consider results that can be obtained after collecting an integrated luminosity of 
200 ft) -1 , corresponding to approximately one year of running at TESLA and to a few 
years of running if parameters proposed for other linear collider options (such as JLC or 
NLC [p9f| ) are used. 

We should stress that one of the highly desirable features of a LC is the ability to tune 
the c.o.m. energy to explore thresholds with precision. In this way, specific signals, e.g. 
from SUSY, can be enhanced from among others, unless the production thresholds are 
too closely degenerate. For instance, we use this property below for our GMSB models 
# 1-2 to isolate the neutralino pair production process for individual study. In addition, 
the energy can be tuned to alter appreciably the Lorentz (flj) factors of the produced 
neutralinos and hence extend the range of NLSP lifetime measurements. Neither of these 
options will be available at the LHC. 

A further advantage to this study of a LC over the LHC is the fact that the effective 
c.o.m. energy is known precisely, up to effects of initial state radiation (ISR) and beam- 
strahlung. The values of <5 p bcam due to beamstrahlung are estimated to be 2.8% and 4.7% 
for the first and second phase of TESLA respectively |JT|. For the processes we study in 



this paper, it is of prime importance to know the energy of the pair-produced neutralinos 
in order to be able to reconstruct the neutralino decay length, as described further in 
Sec. [7]. We wish to stress here the complementary nature of the LC with respect to the 
LHC and we envision the pleasing scenario where the LHC provides a wealth of inter- 
esting data, which is subsequently investigated at the LC with high precision. This may 
indeed be necessary in order to distinguish conclusively between GMSB and other pos- 
sible SUSY realizations (e.g., no-scale SUGRA models) and to measure the fundamental 
parameters with the precision needed to extract striking conclusions concerning physics 
at the messenger and higher scales. Because of this, the present LC detector design pro- 
posals should be flexible enough to benefit at a late stage from LHC new physics data. 
Our study attempts to address this issue and we will further comment on this in Sec. ||. 
Moreover, as we already pointed out in Sec. |], GMSB models often feature very heavy 
strongly interacting sparticles, which could not provide a large signal at the LHC. In this 
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case, the role of the LC in determining the origin of the new physics signal would be even 
more important. 

5 Disentangling a Signal from GMSB at the Linear 
Collider 

In this section, we provide an example of how it might be possible to extract a good 
amount of information about the parameter values of the underlying model from the 
observation of an abundant GMSB signal at the LC and a simple threshold scannning 
technique. Let's assume that Nature has chosen GMSB and that Model # 1 is realized. 
Let's also assume for simplicity that y~F is not too large, so that (most of) the produced 
NLSP's decay within the detector. If this is the case, just a few weeks of LC running 
at some initial c.o.m. energy between 200 and 500 GeV would be enough to recognize 
the presence of an evident GMSB-like scenario. Indeed, a copious number of events with 
two 7's and large missing energy would show up, due to the inclusive characteristics of 
the GMSB signal. Further, we will see in Sec. ^ that in most cases it will be possible to 
show that these photons do not point to the interaction region, and hence are likely to 
come from a delayed neutralino decay, since the SM background is essentially zero. (Of 
course, at least in the case of Model # 1, it is very reasonable that at the moment of 
starting the LC operations clear indications for GMSB would have already come from the 
LHC.) Among the two-photon events, there will be many coming from N±Ni production 
featuring no other particles and, if y/s > 300 GeV, many others including (soft) e + e~ 
pairs from selectron-pair production as well. If the c.o.m. energy is even larger, then 
more complex events, many with hadronic activity, would also appear from, e.g., G\G\ 
production. 

A feeling of the situation can be obtained from inspection of Fig. [5]. In Sec. |], we 
stressed the importance of the ability of a LC of tuning the c.o.m. energy to explore 
thresholds with precision. First, one could vary ^fs in big steps and just inclusively count 
two-photon events to get the rough location of the thresholds for the various SUSY- 
production processes (cfr. thick line labeled "TOT" in Fig. |5|). Then, one could focus on 
the individual thresholds, observe more exclusive characteristics of the signal (for instance, 
£r£r production gives events only, with well-defined lepton energy spectra, since 

.R-sleptons will always decay to Ni, and so forth), and vary y^i in finer steps to get a 
precise value of the sparticle masses involved in the corresponding production process. 
For the case of a GMSB model like Model # 1 where about 10 thresholds are present 
below y/s = 500 GeV, it seems reasonable to assume that a 200 fb" 1 run (less than 1 
year, based on the TESLA expected performance) would allow extraction of the light 
masses with errors at the level of fractions of a GeV. Of course, the fine details depend 
on the slope and the magnitude in the vicinity of the thresholds of the curves for the 
various individual cross-sections as functions of ^fs (cfr. Fig. ||). For instance, in absence 
of important t-channel contributions, one would expect a steeper ~ /3 3 behaviour for 
gaugino-pair (fermion) production compared to ~ (3 for slepton-pair (scalar) production, 
so that gaugino masses could generally be determined with higher precision ||32|| . On the 
other hand, in our case the t-channel contributions are important and, in addition, one 
can always imagine to spend more machine time running close to the "harder" thresholds 
and also use other observables (e.g. distributions) to get additional information on the 
spectrum. 

Our intent here is not to simulate fully such a complex study, but to evaluate what 
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could be the sensitivity in determining the GMSB parameters from the knowledge of the 
light spectrum that could come from a roughly uniform threshold-scanning. Based on 
Model # 1 and a total of 200 fb -1 collected between 200 and 500 GeV c.o.m. energies, 
we estimated the following approximate precisions for the sparticle masses: 



A(m Ri ) ~ 0.2 GeV 
A(m e -J ~ 0.2 GeV 
A(m fl ) ~ 0.8 GeV 



A(m^ 2 ) ~ 0.8 GeV; A(m d J ~ 0.1 GeV; 

A(m e ~J ~ 0.2 GeV; A(m Afl ) ~ 0.8 GeV; (9) 

A(m f2 ) ~ 2.0 GeV; A(m h o) ~ 0.1 GeV. 



The assumption on A(m^o) is based on the fact that many e + e~ — > events would 

be observed at the LC if Model # 1 is realized and many other Higgs events would have 
already seen and studied at the LHCQ 



Also, in Sec. |7.2| , we will see that a measurement of the iV\ mass with a precision at 
the level of a few tenths GeV can be easily achieved by looking at the 7 energy spectrum 
from Ni decays. 

We used a home-made computer program called MinuSUSY ||19|| , interfaced to SUSYFIRE 
and Minuit [[34], to perform fits to SUSY-model basic parameters starting from informa- 
tion on the sparticle spectrumQ The program works both with (m)GMSB and (m)SUGRA 
models and in "global" or "local" mode. The "global" mode is intended to determine 
which class of SUSY models and which approximate values of the basic parameters best 
recover the input spectrum. We used this run mode starting from the light spectrum 
of Model # 1 (cfr. Tab. [I]) and found that indeed there is no mSUGRA model that 
can reasonably fit it. Such a spectrum could be recovered with good precision only by 
releasing one or more of the unification assumptions at the GUT scale. In contrast to 
mSUGRA, the minimal GMSB framework allowed us to single out a successful region of 
the parameter space including Model # 1. Once a rough knowledge of the basic parame- 
ter values is obtained, it is possible to run MinuSUSY in "local" mode around these values 
and get optimised values and errors on them based on the input errors on the sparticle 
masses. Basically, one simulates a large number of possible sets of mass measurements 
using a gaussian distribution for the masses around the central values that one would 
get from the chosen underlying SUSY model. Starting from the errors on the masses for 
Model # 1 quoted in Eq. (|9]), we simulated 100 sets of measurements of the light sparticle 
spectrum. Our results for the 100 subsequent reconstructions of the GMSB parameter 
set performed with MinuSUSY are summarized in Figs. [T7], ITHL for M mess and A, N mess and 



tan p, respectively. Here we did not require A" mess to be an integer and considered it as a 
real variable to perform the fits. (Notice that non-integer values of A mess are possible in 
some non-minimal classes of GMSB models, cfr. e.g. Ref. fl2~3| .) A gaussian (+ constant) 
fit to the distributions gives the results shown in Tab. ||. As for the sign of /1, we found 
that the best fits are obtained for \i > 0, as expected. 

Tab. P indicates that it seems possible to determine A and N mess with a precision of 
about 1 part in 10 3 and tan/3 with 1 part in 10 2 , by just using threshold scanning and 

4 Here, we assume that the theoretical error on determining the lightest Higgs mass from any SUSY- 
model input parameters, currently at the level of at least a few GeV |33| will be reduced by that time by 
more detailed calculations. Similarly, we imagine that the theoretical error on the sparticle masses will 
also be brought at a level comparable to the numbers quoted above. 

5 MimiSUSY does not take higher-order corrections to the sparticle masses into account, but these can 
typically be reabsorbed in a redefinition of the basic model parameters and a shift of the starting values 
needed to generate a given spectrum. For our purpose here, however, the precise values of the SUSY 
parameters are not the point, since we are only interested in evaluating the level of sensitivity one could 
reach by using these techniques. We believe our indications in this respect to be found below are still 
valid without taking fine effects into account. 
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Parameter 


Fitted value 




(161 ± 2) TeV 


A 


(76.01 ± 0.08) TeV 


N 

1 ' mess 


0.9994 ±0.0009 


tan j3 


3.50 ±0.03 



Table 6: Results of fits to the parameters of GMSB Model # 1 starting from a possible 
set of light sparticle masses measurements via threshold scanning technique, as described 
in the text. A 200 fb~ l run at the LC is assumed. 




1500 1525 1550 1575 1600 1625 1650 1675 1700 '725 75400 75600 75800 76000 76200 76400 76600 76800 77000 



M mess [GeV] A [GeV] 

Figure 17: Reconstructed values and fit to M mess (a) and A (b) for Model # 1 as a result 
of 100 possible sparticle spectrum measurements from threshold scanning, as described in 
the text. 

sparticle masses as observables and running for less than 1 year at a LC with ^/s < 500 
GeV. Of course, to achieve such an impressive goal, it is crucial that one can count on the 
high-luminosity, such as the option proposed for TESLA. The only parameter that could 
not be determined at a level of 1% or better is M mess , but this is well understandable 
since the sparticle masses depend only logarithmically on it. Better precision could well 
be reached if other observables (total cross sections, distributions, branching ratios etc.) 
were added to the global fits. On the other hand, it must be said that Model # 1 is a 
particularly "easy" model, in the sense that it yields a light spectrum and the various 
thresholds are well separated. It would be much more difficult if the scenario of Model 
# 2 (for which LC energies well above 500 GeV would be needed to extract the GMSB 
parameters) or Model # 3 (with many sparticles almost degenerate) were realized. 
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(a) 



0.9S 0.985 0.99 0.995 1 

N 



ID 12345 
Entries 100 
Mean 0.9998 

RMS 0.3181 E-02 

xVndf 35.80 / 15 
Constant 10.63 
Mean 0.9994 
Sigma 0.8912E-03 



(b) 



ID 


12345 


Entries 


1 00 


Mean 


3.500 


RMS 


0.1740E-01 


x'/ndf 


18.04 / 52 


Constan 


1.825 


Mean 


3.497 


Sigma 


0.31 16E-01 




Figure 18: As in Fig. [T/], but for N mess (a) and tan/? (b). 



6 Event and Detector Simulation 



To generate GMSB events we used a modified version of SUSYGEN 2.20/03 [35|, where 
the 3-body neutralino decays were added and the corresponding kinematical distributions 
were input numerically, according to the discussion of Sec. |3| and the results obtained with 
Gravi-CompHEP, for our reference Models # 1-3. For each GMSB model, the relevant 
input cards to SUSYGEN were calculated with SUSYFIRE, keeping \/~F (and hence the N 1 
lifetime) as a free parameter, subject to the bounds discussed in Sees. [I] and || The 
generated events were then passed to our detector simulation software, BRAHMS P(J . This 



is a GEANT 3.21 |37j code including material and tracking detectors, as motivated by the 
ECFA/DESY CDR The relevant detector components are simulated as follows. 

The beampipe is taken as a tube of beryllium of radius 1.0 cm and thickness 0.14%Xo, 
where Xo is the radiation length. Five layers of vertex detectors (VXD), each of thickness 
0.12%Vo and point resolution of 3.5 jum are located at radial positions of 1.2, 2.4, 3.6, 4.8 
and 6.0 cm, with respective half ^-lengths of 2.5, 5.0, 7.5, 10.0 and 12.5 cm. 

We include an intermediate tracking chamber (ITC) as material with a total of 0.23%X 
and dimensions 12 cm < r < 30 cm and \z\ < 100 cm, but we do not consider this detector 
for the track fit. In the forward and rear directions, we include a forward tracking detector 
(FTD) made of disks of silicon strip detectors each of thickness 300 /im at z-positions of 
40, 50, 120, 140 and 160 cm, with outer radii of 10, 10, 30, 30, 30 cm respectively and 
inner radii of 2.5, 2.5, 10.1, 11.7, 13.3 cm. All elements have r<p resolution of 25 ^m. 

We use a time projection chamber (TPC) as central tracker with inner active radius 
of 38.6 cm, outer active radius of 162.6 cm and active longitudinal half-length of 250 cm. 
The active volume is filled with gas (which we take to be argon) and provides a maximum 
of 118 hit points along a track, each with point resolution of 160 /im in r<p and 0.1 cm 
in z. The inner wall to the TPC consists of a total of 3%Vo of aluminium. This is an 
important source of conversions which we discuss further in our analysis below. 

For the calorimeter part of our simulator, we use an electromagnetic calorimeter 
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Annular Coverage 


1 cos 9\ < 95 


Barrel r- Dimensions (cm) 


172 < r < 210 


Endcap ^-Dimensions (cm) 


280 < 2 < 330 


Energy Resolution (%) 


10.3/ y/E [GeV] + 0.6 


Spatial Resolution (cm) 


4/yjE [GeV] + 2 


Angular Pointing Resolution (mrad) 


50/ ^E [GeV] 


Time Resolution (ns) 


2/^AE [GeV] + 0.5 



Table 7: Calorimeter parameters used in Monte Carlo smearing. 

(ECAL) and assume simple gaussian smearing with resolutions motivated by those in 
the CDR and given in Tab. |7|. 

In addition to the detector resolutions, there is an additional uncertainty in the position 
of the i.p. due to the beam spot size. In the following, we take the beam spot dimensions 



as given by the TESLA machine design parameters at 500 GeV c.o.m. energy [3l| and 
thus apply gaussian smearing to the production vertex of the neutralinos with <j x = 553 
nm , o y = 5 nm and <j z = 400 am.. 

7 Measuring the NLSP Properties and the 

Fundamental SUSY Breaking Scale at the LC 

In this section, we focus on practical methods to measure the NLSP properties and, in 
particular, its mass and lifetime. As discussed in Sec. in GMSB (or in more general 
LESB) models, the NLSP lifetime can be macroscopic and this opens a very important 
window for inspecting SUSY breaking physics, which is not available in HESB models, 
like mSUGRA. Indeed, as Eq. (181) shows, measuring m NL sp and cr NL sp determines the 



fundamental scale of SUSY breaking yF up to the factor B (cfr. Fig. [14]), that can also 
be measured in principle. We refer to the specific case of neutralino NLSP scenarios and, 
in particular, to the typical Models # 1-3 discussed in Sec. 0. 

7.1 Overview of Experimental Techniques 

In the following, we will first describe several techniques that we propose for performing 
such measurements. Depending on the A^i lifetime, these methods involve (and test) 
different parts of the detector, requiring high and somewhat unusual performances. Hence, 
this study could be an important benchmark in the process of designing a LC detector and 
the related simulation software and should not only be seen as limited to SUSY searches. 

We assume here that at least one SUSY production process (i.e. e + e" — > N\Ni) is 
accessible at the LC, so that two neutralinos (plus possibly some cascade decay products) 
appear for each SUSY pair produced. The neutralinos must then decay to a G, based on 
the discussion of Sec. [3|. 

The topology of the N\ decay is sketched in Fig. ITSI where Fig. |T5|(a) shows the case for 



a purely photonic Vi — > jG decay and Fig. |19|(b) shows the cases where a "charged decay" 
(e.g. Ni — > ji + n~G) has occurred (D = 0) or where the 7 has converted subsequent to a 




Figure 19: Topologies of neutralino (Ni) decays, (a) Pure photonic decay; (b) charged 
decay (where D = 0) or photonic decay plus 7 conversion (D ^ 0) 

photonic decay (D ^ 0). In each event, another N± is present. If the process generating 
the neutralinos is simply NLSP-pair production, then the other N\ moves approximately 
colinearly with the one shown, where the acolinearity angle depends only on the ISR and 
beamstrahlung. This information could be used to constrain the event reconstruction, 
but it requires detailed knowledge of the acolinearity angle distribution. For simplicity, 
we treat each neutralino decay as independent in the following analyses. 

We have seen in Sec. |] that in most cases the Ni decays to a G (which escapes the 
detector) and an observable component. The latter can be a photon or a visible // pair 
coming from a "charged decay". In the first case, the photon can be observed via its 
shower in the ECAL. However, in a real detector there is always material between the 
i.p. and the tracking volumes, so a fraction of the photons coming from SUSY production 
will undergo conversion in the material. Using the tracking detectors to reconstruct the 
resulting e + e~ pairs, it is possible to obtain a very accurate determination of the original 
photon energy and direction. 

In the following, we list the observable final states (for each decayed neutralino) and 
introduce the concepts and methods that we will use later for Aq lifetime measurements. 

a) Photon. We use this final state in the calorimeter pointing technique (see Sec. |7.5|) for 
laboratory decay lengths L = /57cr NLgP between approximately 5 cm and 200 cm. 
Calorimeter timing provides additional information for decay lengths from about 
20 cm to about 120 cm (see Sec. |7^| ). Also, the statistical method of Sec. [7TV| 
will be based on counting events with photonic final states. Furthermore, photon 
conversions in detector material can also be used to measure L in the short range, 
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however we will show below that greater precision can generally be obtained by 
using "charged decays" only, after eliminating conversions by using appropriate 
experimental cuts. 

b) e + e~ pairs. "Charged decays" to e + e~ occur relatively abundantly with BR's of order 
1-2% percent (cfr Sec. ^) and we will use them to measure values of L < 20 cm. 
In addition, this final state also occurs at the level of a few percent when photons 



convert in detector material. We will show in Sec. 7^ how to differentiate between 
conversions and "charged decays". For the latter case, the reconstructed e + e" vertex 
corresponds to the decay vertex of the Ni. For the former case, an extrapolation is 
required in order to obtain the neutralino decay point, see Fig. |19|(b). This procedure 
is discussed in detail in Sec. 01 below. 



c) or h + h~ pairs. (Here h is a charged "stable" hadron.) To improve statistics 

slightly, we will use these final states together with b) in the tracking methods and 
consider an inclusive general "two-track" topology. Note that in these cases, the 
events are always a result of N\ "charged decays" and the relevant BR's are typically 
at the level of several percent (cfr Sec. 

In order to reconstruct the decay parameters from observation of e + e" — > N\Ni events, 
we use the formula 

cos0= 7r~Er> 10 

where E is the nominal beam energy corrected for average losses due to ISR and beam- 
strahlung, and po = (Eq — m 2 ^ )a. This formula allows the angle to be determined 
directly from the measured photon energy. The explicit occurrence of rrift in this equa- 
tion emphasises the necessity of a good neutralino mass measurement. 

For the cases where the photon converts to an e + e~ pair, or on the case of a "charged 
decay", the line of the photon flight is determined from the reconstructed vertex and 
momentum of the pair. This line combined with the value of obtained from the energy 
measurement gives an unambiguous value for A [cfr. Fig. |T9|(b)]. 

In the calorimeter pointing method (discussed in detail in Sec. [7.5| below), the angle ip 
and the distance R are determined directly from the calorimeter shower reconstruction. 
In this way, the decay length A is determined on an event by event basis. 

In the calorimeter timing method (discussed further in Sec. |7]6| below), the time mea- 
surement gives the quantity D + A — R, the shower position reconstruction gives the value 
of R and the energy measurement gives the value of (p. Closure of the triangle allows a 
solution for A to be obtained in most cases (up to a quadratic ambiguity which is resolved 
by the requirement that the decay should take place within the dimensions of the ECAL). 



7.2 Measurement by E 1 Spectrum End-Point 

As we have seen, a good neutralino mass measurement is a central requirement to most of 
the lifetime measurement techniques discussed below. Further, precise knowledge of 
is essential to extract the parameter \fF from crft (cfr Sec. |3|). Here we discuss how this 
measurement could be made at a LC. 

One way to measure is by determining the end points of the photon energy 
spectrum from the N± decay. When many SUSY production channels are open in a 
neutralino NLSP scenario, one gets photons from Ni — > G7 decays with a complicated 
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Figure 20: Photon spectrum resulting from 200 GeV neutmlino pair production with 200 
fb^ 1 at ^fs = 500 GeV. ISR, beamstrahlung and detector effects have been included. Also 
shown is the result of the fit described in the text. 



energy spectrum. If the process e + e" — > NiNi is the only one allowed by kinematics, 
then before radiative corrections the lower and upper ends of the E 1 spectrum are always 
given by 



E, 



min,max 



4m 



'111 



In general, while the most energetic photons will always come from those neutralinos that 
are directly pair-produced, the lower end of the spectrum will be degraded by the presence 
of softer photons coming from other SUSY processes, in addition to the SM background 
(cfr. Sec. |7.8| ). For this reason, we concentrate here on the upper end of the spectrum 
to extract the N\ mass. The spectrum which would be obtained after detector effects for 
200 GeV neutralino pair production at = 500 GeV is shown in Fig. |2(] for 200 fb _1 
integrated luminosity (corresponding to a run of less than 1 "year" of 10 7 s). Here we 
have simulated the SUSY signal that one would detect if Model # 2 was realized. In this 
case, we have seen in Sec. ^| that a(NiNi) = 42.3 fb and all the other SUSY-production 
processes would be below threshold at yfs = 500 GeV. 

In order to extract the functional form of the high edge of the spectrum after ISR, 
beamstrahlung and detector effects, a much larger number of Monte Carlo events was 
used to obtain a fit function using the known Ni mass as an input. The functional form 
includes two exponentials to allow for ISR and beamstrahlung together with a cumulative 
normal distribution, Freq, to account for the calorimeter resolution integrated over the 
sharp edge of the spectrum. The function thus obtained was 
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Figure 21: 3D- and end-view of a representative two-track (e + e + 7 + ft) event, fully 
simulated in the proposed LC detector. 



1 dn (E-E^\ (E-E™ x \ 

__ = 0.788- 0.716 exp^^J- 0.0722 exp^J-J 

lg7 7 j , (12) 

where no and E™ &yi are now free parameters. This functional form was then used to fit 
to the 200 fb _1 worth of simulated data as shown in Fig. ^0] to give the result _E™ ax = 



(200.21±0.21) GeV. To extract the neutralino mass from this measurement we use Eq. (|TTf) 
to obtain = (199.7 ± 0.3) GeV. In principle, the N\ mass can also be obtained to 
high precision at a LC by scanning over the threshold region and previous studies [j3~Cfl 
suggest that a precision of order 100 MeV could be obtained from such scans (cfr. also 
Sec. [Sp. However, our point here is that even from an early run with a few months worth 
of data collected at a nominal, fixed c.o.m. energy, a neutralino mass measurement with 
precision at the level of 2% would be possible. 



7.3 Measuring the Neutralino NLSP Decay Length Using 2D 
Projective Tracking 

In this section, we concentrate on determining the Ni average decay length L = f3^cr from 
the distribution of reconstructed two-track events. To this purpose, we require events with 
at least one photon in the ECAL (coming from one of the two neutralinos produced) and 
in addition two charged tracks which can be reconstructed to form a vertex (coming from 
the other neutralino decaying through "charged channels" or to a photon that converts). 

In Fig. |2T], we show a representative two-track event among those we generated and 
fully simulated in the proposed LC detector (3D- view on the left, end-view on the right). 
This particular N\N\ event features one displaced photon and a e + e~ pair coming from a 
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Figure 22: Results of a fit to the projected lifetime as described in the text (a) for Model 
# 1 and (b) for Model # 2. tn both cases the error bars correspond to 200 fb^ 1 

N\ "charged decay". The non-zero impact parameters of all particles are clearly visible; 
ECAL showers and tracks are shown (the invisible 7 path is also indicated). Only the 
vertex detectors and central trackers are displayed. 

The detector track hits are provided by the BRAHMS output and these were then formed 
into tracks and vertices using a home-made reconstruction algorithm. The photon con- 
version algorithms and multiple scattering effects are internal to BRAHMS and so we have 
thus implicitly taken full account of the detector material present. However, no spe- 
cial provision was made for multiple scattering or for pattern recognition effects at the 
reconstruction stage. 

We concentrate first on the case of very short Ni lifetimes, less than a few mm, where 
all the Ni decays take place within the beampipe. In this case, any reconstructed e + e~ 
pairs will be due to "charged decays" only and so there will be no confusion arising from 
conversions. For this region, we must be aware of the beamspot size, which has an rms 
spread in z typically of 400 /im, meaning that a three-dimensional decay vertex is no longer 
useful for lifetime measurements. Instead we must project the decay vertex o nto the x y 
plane and determine the decay length from the resulting distributions of r = yAl + A^. 

In the following we adopt a conservative approach of using the full GEANT Monte Carlo 
to generate event samples for a range of true (input to SUSYGEN) N\ decay lengths and then 
fit the resulting projected lifetime distributions to a simple exponential, plus a constant to 
allow for long tails. The use of this simple fitting function is conservative in that, once a 
specific GMSB model is chosen, the projected decay distribution could be chosen exactly 
and hence fit the data more accurately. Only decay length measurements greater than 
10 /im are used in the fits, so as to eliminate any residual SM backgrounds with tracks 
originating from the i.p. 

For Model # 1, a set of points were generated using 37,600 (corresponding to 200 fb _1 , 
with <r(A r 1 A r 1 ) = 188 fb at y/s = 270 GeV) fully simulated events for each point. Note 
that only the "charged decays" are observed here, because any conversions will take place 
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outside the beampipe, and so the distributions obtained from this sample can be applied 
directly to any model which has a very short-lived NLSP. 

The results are shown in Fig. ^(a) where the outcome of the fits to the projected 
lifetime is plotted against the true value of L. The vertical error bars shown are the 
parabolic errors obtained from the fit and so correspond to what could be extracted from 
a run of 200 fb -1 . The lifetime is read-off from the fit straight line in order to give the 
true decay length directly; note that any effects due to ISR or vertexing systematics are 
included automatically by such a procedure because all these effects are present in the 
Monte Carlo. This of course also applies to the relativistic factor (f3j). Reading from the 
plot, the error bars at the 50 /im point correspond to a 10% error in the decay length 
measurement whereas the error bars at 500 /im correspond to a 4% error. It can be 
seen from the straight-line fit to the points that decay lengths down to 30 /im are also 
well- measured for Model # 1. 

Performing the same fits to the projected distributions according to Model # 2 with 
8,460 events (corresponding to 200 fb~\ with cr(JViiVi) = 42.3 fb at ^/s = 500 GeV), 
gives an error of 13% for a 50 /im decay length and 10% for a 500 /im decay length. The 
results are shown in Fig. |22|(b). It should be noted that, although Model # 2 provides less 
statistics due to a lower neutralino production cross-section, the possibility of Ni decay 
to an on-shell Z° boson with subsequent decay to lepton pairs with significant opening 
angle somewhat compensates for the loss in statistics. 

Fig. 22 also shows that lifetimes shorter than 30 /im could also be measured in principle 
using this technique. However it should be remembered that, while the intrinsic beamspot 
size may be of order 5 nm in y and 500 nm in x, the actual position of the interation 
point may have to be determined on a pulse-by-pulse basis which could lead to additional 
transverse uncertainties. It should also be remembered that we are fitting to track pairs 
with essentially zero opening angle, which means that simple assumptions about vertex 
resolution must be avoided. So in the following we remain conservative in our claim that 
lifetimes as short as 30 /im could be measured. 

Notice that for Model # 1, the minimum value of y/~F allowed by theory is about 110 
TeV (cfr. Sec. ^|), corresponding to a N\ lifetime of about 200 /im and to a decay length 
L ~ 180 /im when running at ^/s = 270 GeV, while for Model # 2, v^min ^ 212 TeV 
=^> cr- 1 ™ ~ 70 /im ==>- L^ =500 GeV ~ 53 /im. So, our result is that at least for these 

two particular cases, this method more than covers the lower end of the range for the N% 
lifetime. On the other hand, we know that neutralino NLSP models exist where ct^ is 
as short as 5 to 10 /im (cfr. Fig. [15] in Sec. |3|). For these cases, an alternative approach 
would be to increase the c.o.m. energy of the machine, hence increasing the boost of the 
neutralinos, which so far have been close to threshold. In this way, shorter decay lengths 
factors are readily accessible and would realistically be the preferred approach to measure 
lifetimes down to 10 /tm (notice that {f3 r ))f }l can never exceed 2.3 for ^Js < 500 GeV). 
On the other hand, we checked that GMSB models with a very short N\ lifetime tend to 
feature a greater degree of degeneracy among the lightest states, so that running at higher 
energies would also often imply opening other SUSY production channels (i?-slepton pair 
production). As a consequence, it is important to address the problem of being able to 
extract a precise lifetime measurement in the presence of a complex SUSY signal, using 
appropriate selections. 

As an example, we study Model # 3 at ^fs = 470 GeV with 50,400 total generated 
events before cuts, corresponding to 200 fb -1 , since a(NiNi) = 124.6 fb, a{e R e R ) = 64.3 
fb, cr(fifi) = 31.7 fb and a^fiRjiR) = 31.4 fb. The events were then selected at the 
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Figure 23: (a) Reconstructed vertex projective radial distances, r, for Model # 3 with 
L = 10 fim together with a Gt to an exponential plus constant, (b) Shows the results of 
the corresponding fits for a range of L values. 

generator level by the following cuts. The events had to contain either exactly two photons 
with |cos#| < 0.95 or exactly one photon with |cos#| < 0.95 and one charged pair of 
tracks each with |cos#| < 0.99. In addition to these cuts, an event which contained any 
additional particle with |cos#| < 0.99 was rejected. Afer these cuts the sample contained 
21,108 NiNi events together with a SUSY "background" of 1,066 selectron pairs, 463 
smuon pairs and 336 stau pairs. These events were then passed to the BRAHMS detector 
simulation followed by the 2D projective fit procedure. 

Fig. |23|(a) shows the distribution of reconstructed vertex projected radii. The remain- 
ing SUSY background consists of only 3 stau events, so its effects are negligible in this 
case. The error bars for L — 10 fim correspond to 9.8 /im. We conclude that for a 200 
fb _1 run at a/s = 470 GeV, an upper limit on L of approximately 20 /im could be set 
at the 2a level. The measurement could of course be improved with higher luminosity 
together with running at higher energy to utilize a larger (^7) factor. If we could run 
with 10 3 fb _1 at 800 GeV, then we could expect a gain of a/5 from the statistics together 
with a gain of 2.2 from an improved (/?7). These would allow a cr measurement to a 
precision of approximately 20%. 

7.4 Measuring the Neutralino NLSP Decay Length Using 3D 
Vertexing 

For decay lengths greater than about 500 /im, the three-dimensional vertex information 
is useful because the beamspot size in z is relatively less important. For this region, 
photons converting in detector material must be taken into account because these events 
will appear very similar to the "charged decays". The vertex position of reconstructed 
tracks from these two processes are shown in Fig. |24] for Model # 1 with L = 10 cm. 
For both the "charged decays" and the conversions, a line of flight can be defined 
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Figure 24: Reconstructed vertex positions for Model # 1 with L = 10 cm for (a) "charged 
decays" and (b) conversions of photons from "neutral decays". The radial coordinate of 
the vertex r is plotted against the longitudinal coordinate z. For case (b) the excess of 
conversions in the TPC inner wall and in the inner masks can be seen clearly. 



by the lepton-pair vertex and the direction of the reconstructed pair momentum. In the 
case of conversions, this line of flight reconstructs that of the photon from the neutralino 
decay. The "charged decays" are also included in this approach effectively as photons 
which have converted instantly. Referring to Fig. [l9|(b), cos0 is now determined by the 
energy of the reconstructed pair and we have 

A = R 8 ^ (13) 
sin <p 

In this way a A measurement can be obtained for conversions as well as for "charged 
decays" . However, because we do not wish to make assumptions on the two-track resolu- 
tion of the central tracker (and hence the reconstruction efficiency for conversions in the 
outer detector), we prefer here to use only "charged decays" for lifetime determination. 



However, as shown in Fig. 25|(a), we can use the variable p = X/R as a discriminant. In 
the following we require p > 0.5 to improve the "charged decay" purity. 

The lepton-pair invariant mass distributions for "charged decays" and conversions are 
shown in Fig. [251(b) for Model # 1 with L = 10 cm. Both peak at very low values, of 
order 10 MeV, for the detector resolutions employed here. However, the conversions are 
more sharply peaked than the "charged decays", so cutting on 20 MeV < M pair < 10 
GeV improves the purity of the "charged decay" sample. In addition this cut is desirable 
because it selects events with a larger opening angle between the daughter tracks, which 
results in an intrinsically more accurate vertex position measurement. The higher mass 
cut is to eliminate any residual SM background from leptonic IV-pair events. 

In addition to these cuts we make an additional simple geometrical cut of r < 30 cm 
to remove the large number of conversions in the inner TPC wall. After performing these 
three cuts, the A distribution is primarily from the "charged decays" and this is illustrated 
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Figure 25: These figures are for Model # 1 with L = 10 cm with statistics corresponding 
to 200 fb^ 1 . The unshaded histograms are the conversions and the shaded histograms are 
the "charged" decays, (a) Shows the p parameter distribution and (b) shows the invariant 
mass of the reconstructed track pair. Note the logarithmic scales. 



in Fig. |26], where the fit to Model # 1 with L = 10 cm is shown. It is clear from this 
figure that conversions have been reduced to a negligible level. 

This procedure was repeated for a range of L values from 500 /im up to 15 cm. The 
results are shown in Fig. |27] where each point represents the result of the exponential fit, 
with corresponding errors for samples consisting of 37,600 generated events (200 fb _1 ) 
before cuts. The Monte Carlo includes the effects of ISR. 



Also shown in Fig. ^7] is a straight line fit to the data in (a) which is extrapolated 
to (b). Using this fit we obtain a relative measurement statistical error of 5% for L = 1 
mm and a relative error of 4% for L = 10 cm. Our cuts are not optimised for L greater 
than about 10 cm. In this region geometrical detector effects and variations in track 
reconstruction efficiency become important, as evidenced by the deviation from a straight 
line in (b). Rather than correct for these effects in detail here, we note that this region 
of L is very well covered by the calorimeter pointing method discussed below. 

We conclude that it is possible to make a good measurement (to 10% or better) of nom- 
inal decay lengths ranging from 30 /im up to 10 cm under very conservative assumptions, 
using the tracking detectors. 



7.5 Measuring the NLSP Decay Length Using Calorimeter 
Pointing 

A finely segmented ECAL allows for the possibility of detecting photon impact parameters 
with respect to the i.p. The direction of the photon is reconstructed by fitting to the 
distribution of energy deposits among the individual ECAL cells which make up the 
electromagnetic shower. The direction finding is improved significantly by the use of 
presamplers to provide a precise point along the photon direction. A typical presampler 
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Figure 26: Distribution of decay length, X, for Model # 1 with L = 10 cm after 200 fb . 
The 3D vertex reconstruction is used and the distribution is fit to an exponential plus 
constant. The contributions from conversions remaining after cuts are included in the 
unhatched histogram, but are also superposed separately as a hatched histogram. 




Figure 27: Results of fits to the 3D decay lengths as a function of the true value of L 
using the tracking method for Model # 1 . Part (a ) for 500 fim < L < 1 cm and (b ) is for 
the results for 1 cm < L < 15 cm. The error bars correspond to 200 fb' 1 . The fit lines 
are described in the text. 
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Figure 28: Results of fits to decay lengths derived from pointing information. The data 
points and error bars correspond (a) to a nominal 37,600 neutralino pairs for Model # 1 
and (b) 8,460 events for Model # 2 (corresponding to 200 fb~ x for each model) . The 
curves correspond to quadratic polynomial fits to the data for input neutralino masses of 
±100 MeV from the true mass. 



detector consists of lead/scintillating fiber structures which both initiate the shower and 
measure the point at which the shower begins to a precision of a few tens of /xm. The 
use of more than one presampler can provide two points along the photon direction and 
further improve the direction finding. 

In the following, we assume that the angular resolution of the ECAL is that given 
in Tab. |7|. We also need to take account of the spatial and energy resolutions as listed 
there. Referring to Fig. 0(a), the angle and the distance R are determined directly 
from the calorimeter shower reconstruction. The angle (f) is calculated from the energy of 
the shower, together with the measured value of the neutralino mass, as described above. 

The effect of the mass measurement error was estimated as follows. Each Monte Carlo 
event was reconstructed using the nominal 100 GeV mass and also using assumed mea- 
sured masses of 100.1 GeV and 99.9 GeV. We showed in Sec. [T2]that a mass measurement 
error of 0.3 GeV could be obtained in a first run of the detector using the upper end of 
the photon energy spectrum. It would be natural to assume, given any such discovery 
of neutralino production, that sufficient luminosity would be made available to perform 
detailed threshold scans and reduce the mass uncertainty to the level of 0.1 GeV (cfr. 
Sec. |). 



The results are shown in Fig. |28|. The Monte Carlo includes effects of ISR and each 
point was obtained using 37,600 generated events before cuts for Model # 1 and 8,460 
events for Model # 2, corresponding to 200 fb _1 for each model running at ^Js = 270 
and 500 GeV, respectively. At each point, the reconstructed decay lengths A were fit 
to an exponential for 5 < A < 200 cm and the points show the fit decay length with 
corresponding error. The A > 5 cm cut, corresponding to about 2<r for the resolution of 
this method, serves to reduce any SM background to a negligible level. 
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The resulting calibration curve is well approximated by a quadratic polynomial and 
the corresponding fit curves obtained assuming an error on the neutralino mass of ±100 
MeV are shown together with the points obtained using the correct input masses of 100 
GeV for Model # 1 and 200 GeV for Model # 2. In this way the systematic error due to 
the mass uncertainty is shown to be very small for L less than 200 cm for both models. 
If we examine the lower statistics Model 2 in the region of L = 200 cm, we find a 
statistical error of 6% and in the region of 10 cm we find a statistical error of 1%. Our 
conclusion is that the calorimeter pointing method works very well for 5 cm < L < 2 m. 

7.6 Measuring the NLSP Decay Length Using Calorimeter 
Timing 

Calorimeter timing information is highly desirable in any detector at a LC in order to 
reject cosmic rays and many beam-related backgrounds as well as for its use in a trigger. 
In this section, we investigate a further use of sub-nanosecond timing for neutralino decay 
length measurement. 

For the purposes of this study, we use the energy-dependent timing resolution given 
in Tab. [7]. Referring to Fig. |19|, a calorimeter timing measurement gives the quantity 
D + A — R. This timing information is then combined with the position and energy 
measurements, all smeared according to the resolutions given in Tab. [7|. 

The timing resolution that we use is motivated by what could realistically be achieved 
for a large calorimeter. This is of the order of 0.5 ns, which naively implies an intrinsic 
photon impact parameter resolution of order 15 cm. Data sets for a series of decay lengths 
were simulated using the full Monte Carlo and the reconstructed decay length distributions 
were fit to a simple exponential between the lengths of 30 cm and 120 cm. The lower 
limit is given by the timing resolution and the upper limit by the geometrical acceptance 
of the detector. The results of the fits are shown in Fig. |29|. The data points correspond 
to 200 fb _1 (37,600 and 8,460 generated neutralino pairs at the usual c.o.m. energies for 
Models # 1 and # 2 respectively, before any cuts are applied). The decay lengths are 
determined on an event by event basis assuming the exact neutralino mass. The effect 
of a ±100 MeV shift in the input raff is indicated by the curves, which were obtained 
by fitting the resulting data points (not shown) to a quadratic polynomial. The results 
confirm that the effect of neutralino mass uncertainty on a lifetime measurement is small. 
Values of L below 10 cm are poorly determined due to the intrinsic timing resolution and 
large values of L are limited by statistics; the worst case is for Model # 2, where statistics 
get poor for L > 120 cm, so we take this as the upper limit for this method. 

This study shows that timing information alone can not improve on what could be 
achieved using calorimeter pointing and tracking methods. The timing information could 
of course be combined with the tracking and pointing methods to achieve an improved 
accuracy, but we do not consider this further here. We stress, however, that timing should 
not be neglected in the overall detector design because it may be needed to reduce cosmic 
ray backgrounds and beam related backgrounds as well as being of use in the trigger for 
GMSB events. We will further comment on this in the next sections. It is possible that 
the NLSP is so long-lived that the photon decay would appear in the ECAL several bunch 
crossings later than expected and indeed not in time with any bunch crossing. To detect 
such events good timing information would clearly be required. 
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Figure 29: Results of fits to decay lengths derived from timing information. The data 
points and error bars correspond (a) to a nominal 37,600 neutralino pairs for Model # 1 
and (b) 8,460 events for Model # 2 (corresponding to 200 fb -1 for each model) . The 
curves correspond to quadratic polynomial fits to the data for input neutralino masses of 
±100 MeV from the true mass. 



7.7 Measuring the NLSP Decay Length Using Statistics 

A measurement of the NLSP lifetime can also be made using a simple counting technique, 
because the probability that, e.g., a photon from the N\ — ► 7G decay is observed in the 
ECAL is a function of the N\ lifetime. As a consequence, the ratio of two-photon to 
one-photon events observed after SUSY production is a function of the ct^ . 

As anticipated in Sec. |] for the spherical detector case, this function is in principle 
a simple combination of exponentials convoluted with the effects of detector geometry, 
cuts designed to eliminate SM backgrounds and calorimeter performance. We determine 
the functional dependence of the ratio using Monte Carlo techniques. Neutralinos were 
generated from pair production using our modified version of SUSYGEN interfaced to a 
modified version of CIRCE |38j to include full effects of ISR and beamstrahlung. Any 
photons coming from N± — > 7G decays with energy loosely included in the range given 
by Eq. 



11 



and which originated within the detector tracking volume were extrapolated 
and required to hit within the acceptance of the ECAL. Any event with more than two 
photons or less than one photon was rejected at this stage. We also imposed cuts on the 
missing energy, where appropriate, based on the expectations for the signal. 

The position of the photons and their energies and angles of pointing were smeared 
according to the parameters in Tab. |7]. The resulting photon impact parameters at the i.p. 
were then calculated and at least one photon in the event was required to have an impact 
parameter greater than 30 cm. If one could believe in gaussian statistics to this level, this 
would correspond to a 10cr cut and is hence designed to reduce any SM background to a 
negligible level (we will further comment on this later in this section). 

We define a one-neutralino event as containing one large impact parameter photon and 
nothing else, and the total number of these events is . We define a two-neutralino 
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Figure 30: Plot (a) shows the statistical ratio functions for Models # 1 and # 2, at 
\fs = 270 and 500 GeV respectively, as discussed in the text. Plot (b) shows the decay 
length reach for the same models as a function of integrated luminosity; the decay lengths 
which can be measured to an accuracy of 10% are bounded by the curves. 

event as one where there is at least one large impact parameter photon together with 
anything else visible in the detector, and the number of these events is n 2 ^ . Note that 
n 2JVi does n °t include two-photon events only, but also events where both neutralinos 
decay visibly inside the detector; one of them through N-y —>■ 7G and the other through 
any channel. Our choice is designed to maximise statistics, but we note that simple 
experimental cuts based on the event topology could easily remove any event which is 
not purely photonic, should that yield a cleaner analysis. We checked that this would not 
result in an important loss of statistics (up to about 15% on n 2 ^ in the worst cases with 
a heavy N\), so that the essential conclusions of our analysis would remain unchanged. 
As for n 1] y , instead, we only include photonic decays because we rely on the presence of 
a non-pointing photon to eliminate the irreducible radiative SM background, as discussed 
below. 

The ratio R<j = — ^ was then determined as a function of the neutralino decay length 

for Models # 1 and # 2 by generating a nominal 10 7 events running at ^/s = 270 and 
500 GeV, respectively, to start with. The resulting functions are shown in Fig. [301(a). 

In order to extract the precision AL on the measured decay length, we take the error 
Ai?°- 7 on the measured ratio for the 10 7 events to be: 

iVi 



and then we obtain 



AL 1 /10 7 \^ . „ (dR Nl 



(14) 



-1 



A< (15) 



L U to J ^ I dL 
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Figure 31: Figure (a) shows the optimal c.o.m energy at which to run the LC as a function 
of ct for Model # 2. Figure (b) shows the relative precision of ct obtained after a run of 
200 fb -1 at the corresponding optimised c.o.m energy. 



where n tot is the total number of expected events. Note that n tot = JCa(NiN\) dt is not 
equal to the sum of and n 2 fj 1 because some events will not be counted, due to both 
neutralinos decaying outside the detector acceptance and/or to non-selected final states. 

The resulting curves giving the decay length reach are shown in Fig. |30|(b). Assuming 
an integrated luminosity of 200 fb -1 , the plots show that laboratory decay lengths are 
well measured by this technique down to approximately 50 cm for both models and that 
the upper limit ranges from approximately 20 m for Model # 2 up to 40 m for Model 
# 1, corresponding to cr values of 27 m and 44 m, respectively. It should be noted that 
the sensitivity of the statistical technique is dependent on the c.o.m energy. The statistics 
are determined by the quantity L = fi'ycT, whereas the parameter of physical interest is 
cr. We have seen that if we take Model # 2 with c.o.m. energy 500 GeV and allow a 
maximum measurement error (4^ ) =0.1 then, as indicated above, we find that we 

\ J ' / max 

can measure L up to 20 m, corresponding to cr max = 27 m. If we now decrease the c.o.m 
energy, then we lose in cross section, but gain in intrinsic sensitivity. In this way, an 
optimal c.o.m energy can be found for each value of cr and this optimal energy is plotted 
against cr for Model # 2 in Fig. [H](a). The corresponding optimised precision A ^ is 
shown in Fig. |3l|(b). Unfortunately, it turns out that, at least for Model # 2, any overall 
gain after optimisation is small. The larger value of cr that can be measured with a 
precision of 10% or better is again about 27 m, and the corresponding "optimised" c.o.m. 
energy is about 463 GeV. While this shows that optimisation gains tend to be small, it 
also tells us that the method is not critically sensitive to the c.o.m. energy, even though 
we are running in a threshold region. 

To evaluate the level of background reduction we need for a meaningful analysis, let 
us consider Model # 2 with cr = 27 m and optimise the c.o.m. energy to 463 GeV. 
At this energy, the corresponding laboratory decay length is L = 14.6 m and the NiNi 
production cross section is 25.5 fb. Using 200 fb -1 , we obtain the expected values of 
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n iJVi = ^25 and n 2 fj = 178. The irreducible backgrounds would need to be subracted 
from the raw counts involving numbers of this order of magnitude, so we need to check 
that after a suitable set of cuts the number of background events left is small compared to 
the statistical errors of about ~ 29 for one neutralino events and vl78 ~ 13 for two 

neutralino events. Experience at LEP (see for example pSfl), has shown that backgrounds 
from cosmic rays and detector noise (sparks and radioactive decays) can be reduced to 
negligible levels compared to the physics backgrounds by requiring the event timing to be 
consistent with a beam- crossing. Leaving aside the minor case of events counted in n 27 y 
not including 2 photons (they are easily treatable anyway), the SM physics backgrounds 
are expected to be dominated mainly by 1/1/7(7) events and, to a lesser extent, by radiative 
bhabha's. It thus remains to check that the cuts applied in the counting procedure reduce 
the expected number of 1/1/7(7) events to a "few" events. For this purpose, we performed 
a quick evaluation of the cross sections for these processes using CompHEP^ and requiring 
the following for each photon: a) |cos# 7 | < 0.95; b) E 1 included in the range flTTI). We 
then found that: 

a cvT{e + e" -> ViVfi) ~ 810 fb at y/l = 270 GeV; 

i 

570 fb at Vs = 500 GeV; 
<TcuT(e + e~ — > ViVijj) ~ 31 fb at y/s = 270 GeV; 

i 

7 fb at y/s~ = 500 GeV. 

On top of this, one has to take the 30 cm cut on the non-pointing photon impact 
parameter into account. Of course, it is not reasonable to trust the tails of gaussian 
statistics so much as to consider this as a real 10cr cut. Nevertheless, we note that if 
we were to "conservatively" evaluate its reduction effect and assume it to be equivalent 
to an effective 4cr cut, the one photon background in 200 fb -1 would still amount to 10 
(7) events only for Model # 1 (2), while the two photon background would basically 
disappear. Hence, we are confident that our set of cuts is strict enough to reduce any 
background to the level needed for our study to be valid. 

While we are aware that detailed studies on non-gaussian tails for calorimeter angu- 
lar pointing would be highly desirable in order to solve more precisely the background 
subtraction problem in such a statistical study, we also note that further background 
reduction is possible before appealing to photon non-pointingness. First, with reference 
to the most dangerous vv^j background, we checked that, after our upper cut on the 
photon energy in Eq. (|TT|) , the remaining contribution from e + e~ — > ^(Z* — > vv) is small 
and more than 70% of the cross section at y/s = 270 GeV comes from pure H^-exchange 
(even more than that at y/s = 500 GeV). On the other hand, jy-exchange graphs tend to 
produce relatively soft photons predominantly along the beam direction, while the N±Ni 
signal gives a flat _E 7 spectrum between the end points ([11]) and does not show an im- 
portant angular structure, due to the isotropy of the Ni decay. Based on this, we found 
that, e.g. for the case of interest for Model # 1 with y/s = 270 GeV, imposing stricter 
cuts E 1 > 40 GeV and |cos# 7 | < 0.9 would result in a factor 2 to 3 gain on the signal to 
background ratio, while the corresponding loss of statistics would be limited to less than 
30%. 



6 At this stage, we did not take into account ISR, beamstrahlung effects or any other effects coming 
from undetected radiation 
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Second, and most important, there is here a good chance of exploiting the beam 
polarisation option at the LC. If, on top of the kinematical cuts mentioned above, one 
could benefit from, say, an 80% polarisation for the electron and a 60% one for the 
positron, one would get an additional 0.3 or so reduction factor on the dominant W- 
exchange contributions to the one photon background, while the NiNi signal would get 
an enhancement factor, since an important part of it comes from i?-selectron exchange in 
the t-channel (we remind that our neutralinos are mostly bino's). If even stronger beam 
polarisation was available, then one would probably be able to reduce the one photon SM 
background to an acceptable level by using optimised kinematical cuts and a looser cut 
on the photon impact parameter such that gaussian statistics could still be trusted. 

On the other hand, we checked that the results we obtained for the upper end of our 
statistical method reach on cr^ is not very sensitive to the cut we impose on the photon 
impact parameter in the 10's of cm range. The lower end of the cr reach does, but we 
know from the previous section that the L = 10 cm-1 m range is well covered by the 
direct calorimeter pointing technique. All these considerations make us confident that 
our statistical study is meaningful and our claim that a measure of cr with a precision 
at the level of 10% is possible for typical GMSB models up to several tens of metres by 
using such a method is safe and possibly even slightly conservative. 

Finally, we would like to comment on the fact that larger values of cr might be 
accessible by considering a statistical analysis based on the ratio R'^ = %jVi/ n oiVi °f one 
(non-pointing) photon events to events where both neutralinos decay outside the detector. 
Of course, in order to count the latter events, it is necessary that other SUSY processes 
are within kinematical reach in addition to N±Nx production, so that events including the 
visible products of the decays to the NLSP are present. Taking Model # 1 as an example, 
one could think of running the LC at y/s = 350 GeV or so, allowing i?-slepton pair 
production followed by £r — > £G (£ = e, //, r) decays with 100% BR. The resulting SUSY 
signal would then be made up mainly of (77)$', e + e"(77)^, fi + f/T (77)^, and r + r~ (77)^ 
events, where one or both of the non-pointing photons in parentheses coming from N\ — > 
7G may or may not originate within the detector, depending on cr^ and the N% boost for 
the specific process. For very large neutralino lifetimes, events including two photons are 
very rare and can be neglected. Note that running at higher c.o.m. energies would increase 
the (^7) factor for AVs coming directly from pair production, resulting in a reduction 
of 

n iN 1 fr° m this source, whereas the other processes would produce softer neutralinos 
with a relatively larger probability of decaying within the detector and contributing to 
n^. In addition, the latter processes would provide all the visible SUSY events to be 
counted as Hence, the dependence of R'^ on cr would now be affected by the 

sparticle spectrum and model details, which renders this option less general and more 
involved than the statistical method based on i?^ = n>2N 1 / n iN 1 - ^ n addition, while the 
SM background to events counted in n 17 y can still be reduced by requiring large impact 
parameters for the single photon, such a drastic procedure is not available to reduce the 
SM contribution to n oi y . In the absence of a detailed mSUGRA-like analysis, this would 
lead to a reduction in the R'^ sensitivity to cr. By quick inspection of Fig. |16| for Model 
# 1 at \/s = 270 GeV, one can get a feeling of the dependence on L. Although 
one should take into account the precise definition of n-yy and the differences arising 
from going to e.g. yfs = 350 GeV from 270 GeV, one can still estimate that neutralino 
lifetimes of the order of 1 km or possibly more might be measurable with some precision 
by using the R'^ statistical method. To further increase statistics and the reach in cr, 
one could even consider going to the highest available c.o.m. energies and include all 
possible processes and cascade decays in the counts. However, this would involve very 
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complex analyses which might be desirable only at a later stage, once it was clear that 
either the value of y/~F to be measured is much larger than the bound suggested by the 
simple cosmology condition rnQ < 1 keV or that only an upper limit on CTft can be set. 
A statistical analysis along these lines, as well as other studies concerning measurements 



of GMSB parameters at the LHC will also appear soon ED 



7.8 Backgrounds 

In this section, we address again the issue of background subtraction and summarize our 
strategy in this respect for the various techniques described above. 

Starting with the very short N\ lifetime and 2D-tracking method case, we note that 
there are no backgrounds from conversions, because we include only those vertices that 
occur well within the beampipe with r < 1 cm (r being the radial coordinate of the 
vertex) in the exponential fit. Also, any tracks originating from the i.p. should not 
contaminate the signal, since we impose a cut r > 10 /zm when fitting. We also require 
the invariant mass of the pair of tracks be either less than 10 GeV or in the (91 ± 5) 
GeV range (to retain the Z° peak in the case of Model # 2) to remove any chance of 
trace backgrounds from leptonic W^-pair decays with a vertex reconstructed at a distance 
greater than 10 /im from the i.p., due to resolution effects. We also require both leptons 
to be of the same flavour, which further reduces this background. As for the Z°{-^ 
up)Z°(— > ££) background, we note that in this case the large 11 invariant mass requires a 
significant opening angle between the tracks, which makes our 10 /im cut on r very severe. 
Any remaining background must then originate from particles with lifetime, such as r's. 
However, r's can only mimic our signal if they decay into 3- (or 5-) prong channels (with 
15% BR), which are also mistakenly reconstructed as 2-track events. The effect should be 
small, but any final detector design should be checked to have sufficient 2-track resolution 
to ensure that it is negligible. We did not address this and the above problems in detail, 
because in our study we can safely require the additional presence of a hard [E 1 > 20 
(45) GeV for Model # 1 (2)] photon coming from the decay of the other neutralino and 
more than 40 (90) GeV of missing energy, together with selecting only those events where 
there is nothing else apart from the photon and the reconstructed track pair. 

As for the short decay length and 3D-tracking method, the only additional feature 
compared with the previous case, is the presence of 7 conversions that may or may not 
be used for the lifetime measurement. We conservatively decided to remove them and we 



describe our procedure for this in detail in Sec. |7^. Note also that our cuts here are on 
the full 3D decay length A rather than on the projection r. 

For intermediate decay lengths (calorimeter pointing/timing methods), the main back- 
ground arises from the 77//^ process. All the photons arising from this channel will point 
to the i.p. and so can be reduced by cutting on the reconstructed decay length A. In 
Sec. [7.5| , we achieved this by starting our exponential fits at A = 5 cm, which corresponds 
to roughly a 2a cut for the pointing accuracy assumed in our study. In Sec. |7.7| , we found 
that Y,i vdlViVi) at y/s = 270 GeV (for Model # 1) is roughly 31 fb after appropriate 
cuts, which means that about 1.5 fb are left after the 2a cut. This implies that the re- 
maining background is less than 1% of the signal and so we are safe in our claim that a 
measurement to 10% or better in the lifetime measurement is possible using these meth- 
ods. It should be noted that we have assumed gaussian angular resolution for photon 
pointing and any effects due non-gaussian tails have been neglected. Such residual effects 
in angular pointing should be included in any future detailed calorimeter studies for the 
LC detector. 
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For large decay lengths using the statistical method, we have to address the problem of 
both 77W and the more severe ^vv backgrounds. Any backgrounds due to SM processes 
will be independent of the neutralino decay length and so the SM contribution can be 
estimated from Monte Carlo and be subtracted from the measured values of n x ft (and 
n 2 ft) i n our curves in Fig. |3U], provided the subtractions are relatively small. To this pur- 
pose, in addition to cuts on the photon and missing energies designed to fit the spectrum 
of the signal, we require a very strict 30 cm cut on the non-pointing 7 impact parameter 
to reduce the 7W background to a negligible level. This is again somewhat dependent on 
neglecting non-gaussian tails in the calorimeter resolution. However, in Sec. |7]7], we also 
discussed in detail possible alternative strategies involving kinematical cut optimisation 
and beam polarisation. 

There remains the question of backgrounds from cosmic rays and beam-related back- 
ground. As we have mentioned in Sec. [7.7| , experience from LEP has shown that cosmics 
can be rejected provided one can apply good timing information and require no hadronic 
energy in the event. In addition, the background from cosmics (and beam- related back- 
grounds or calorimeter noise) can be estimated using random dedicated triggers and sub- 
tracted if necessary. Additional veto walls could be constructed around the detector to 
kill any remaining muonic background, if the above measures proved insufficient. 



7.9 Summary 

We have shown in this section that a general purpose detector operating at a LC with 
c.o.m. energies in the 200-500 GeV range could provide a good N\ lifetime measurement, 
to 10% or better, over 6 orders of magnitude, from tens of fim to tens of m for a few 
representative scenarios. We have based this on a 200 ftr 1 run, corresponding, e.g., to 
approximately 1 year at the high luminosity proposed for TESLA and a few years at 
JLC/NLC. 

The techniques we have used in the study are summarised in Fig. ^2[ as a function of 
the Ni average decay length. It can be seen that while a wide L range (e.g. 30 /zm-40 m 
for Model # 1) can well be covered by using individual techniques, there is also significant 
overlap for many of the intermediate regions. This implies important redundancy in the 
lifetime measurement and the opportunity to combine results to achieve an even greater 
precision. 

For many cases (e.g. Model # 1), this is more than enough to cover the lower end of 
the theoretically allowed range for crft and the upper end as well, if naive cosmological 
constraints are imposed. On the other hand, we recognise that scenarios with very short 
lifetimes at the level of a few fim are possible in the context of GMSB and in this case one 
can either set an upper limit or push for higher LC energies to maximise the relativistic 
factor (flj) ft (facing however more complex analyses in the presence of competing SUSY 
signals). 

We have also considered scenarios where the gravitino mass is heavier than ~ 1 keV 
and consequently the neutralino has a very long lifetime > 10's of metres. For these cases, 
we have proposed a way to slightly extend the cr reach by running the LC at a c.o.m. 
energy optimised with respect to the relativistic boost factor and a simple statistical 
analysis, but we have seen that a solid improvement can only come from increasing the 
available integrated luminosity. Interesting prospects for a measurement of very large 
values of or ft up to ~ 1 km also exist, based on a more complex statistical analysis in 
a SUGRA-like scenario where most produced NLSP's appear stable and are invisible, 
but there is still a non-negligible number of neutralinos undergoing a decay within the 
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Figure 32: Summary of the various techniques we have proposed to use at a LC for a Ni 
lifetime measurement with a precision at the level of 10% or better. 



detector. 

Referring to Eq. @, we note that a 10% error in cr corresponds to a 3% error in 
y/F. This is of the same order of magnitude as the theoretical uncertainty on the fac- 
tor B introduced in Sec. |3|. In comparison, the contributing error from the neutralino 



mass measurement using the threshold-scanning technique or the method of Sec. |7.2| is 
negligible. 

Hence we conclude that, for the models considered and under conservative assump- 
tions, it is possible to determine y/~F with a precision of approximately 5% by only per- 
forming N\ lifetime and mass measurements in the context of GMSB with neutralino 
NLSP. Less model dependent results can be obtained by adding information on the Ni 
physical composition from other observables, such as Ni decay BR's, cross sections etc. 



8 Conclusions 

After introducing the GMSB framework and discussing the region of the parameter space 
of interest for LC searches, we focused on the case of a neutralino NLSP and demonstrated 
how measurements made at the LC can provide information on the detailed structure of 
the theory at both the electroweak and the very high energy scales. 

We have shown how a study of the SUSY particle mass spectrum, measured e.g. via 
threshold-scanning techniques, can allow the determination of the fundamental GMSB 
parameters with high precision. In an explicit example where many SUSY thresholds can 
be explored running at y/s < 500 GeV, we found that accuracies at the level of 0.1% for 
A and A mess , 1% for tan/3 and 1-2% for M mess are achievable. 

In particular, we stressed the possibility of performing a measurement of the Ni lifetime 
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at the LC with the aim of extracting information also on the SUSY breaking sector of 
the theory and its fundamental scale y/F. To this purpose, we studied main and rare Ni 
decay channels in detail and set up specific simulation tools based on the TESLA/CDR 
proposals. Using representative GMSB models, we then found that cr^ can be measured 
to 10% or better over a large range, from tens of /im to tens of m, which in many 
cases covers the scenarios allowed from the theory and suggested by simple cosmological 
arguments. (We also sketched a possible way to treat the case where 10's m < cr^ < 1 
km.) We showed that this, together with a Ni mass measurement to ~ 0.1% (whose 
feasibility at the LC we also demonstrated), yields a 5% determination of \fF under 
conservative assumptions, if minimal GMSB-model constraints are used. Better and less 
model dependent results can also come from analyses of other observables such as N\ 
decay BR's, cross sections etc. 

Our results in this respect and the reach in ct^ i depend on details of the detector de- 
sign, which we assumed to be general-purpose oriented. For instance, for intermediate N\ 
lifetimes, we obtained a very high accuracy using calorimeter pointing/timing techniques. 
It should be noted, however, that our calorimeter pointing precision has assumed the 
presence of pre-shower detectors, while simple calorimetry alone might not be sufficient 
for our purpose. 

On the other hand, extreme cases allowed by the GMSB framework with very short 
(few fim) or very long ( > 1 km) Ni lifetime require special detector design. For the very 
short case, one should consider the feasibility of upgrading the interaction region to include 
smaller beampipe radii combined with ultimate vertexing technology. Better performances 
for the very long case can only be achieved by having electromagnetic calorimetry at larger 
distances from the i.p.. This points to either a larger detector or dedicated additional 
devices, e.g. lead scintillator arrays or similar, possibly well separated from the main 
detector (a proposal in this direction for the DO experiment at the Tevatron has already 
been discussed in Ref. If clear SUSY signals are detected at the LHC and/or during 

the first phase of LC operations, then it should be possible to distinguish scenarios with 
a stable N± (e.g. mSUGRA) from GMSB scenarios with a very long-lived N\ NLSP and 
large yF by measuring and studying the sparticle spectrum or other observables. If 
the indications favour a GMSB pattern, then the addition of dedicated devices at large 
distances would be highly desirable for later LC runs. In this case, the position of the 
detector in the experimental hall as well as the dimensions of the hall itself should be 
cleverly designed to allow for such improvements. 

Our study clearly shows that high luminosity at the LC, such as that proposed for 
TESLA is very desirable in order to allow both a better accuracy in extracting the GMSB 
parameters from detailed knowledge of the sparticle spectrum and a good precision mea- 
surement of the NLSP lifetime in (most of) the range suggested by theory. Indeed, we 
checked that any reduction in luminosity would significantly eat away at both ends of 
the cr^ reach, while the situation would improve perceivably if even higher luminosities 
could be available. 

We have pointed out that the clean environment and the flexibility of an e + e" collider 
are ideal to achieve precise measurements in this respect in a variety of GMSB scenarios. 
In particular, we made strong use of the LC ability of running at a variable c.o.m. energy to 
benefit from SUSY thresholds for both measuring the GMSB spectrum and for facilitating 
the N\ lifetime measurement. Also, for the latter measurement, the precise knowledge of 
the N\ production energy was essential and we note that this technique is not available 
at a hadron collider. 
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Finally, we would like to stress that the present study turned out to be an ideal 
benchmark for developing software intended for the TESLA/CDR detector simulation 
that is expected to be broadly used for future analyses of general interest. We have 
indeed seen that all parts of the detector are potentially involved in the measurement of 
the NLSP lifetime and unusual and extreme performance and precision have been often 
required, as e.g. in the analysis based on tracking for the very short cr case. Again, we 
note that the usefulness of the software and algorithms we developed for our purposes 
in GMSB is in many cases extendable to other new physics scenarios, such as general 
LESB models or .R-parity violation. Also, our analysis could hopefully trigger new, non- 
standard ideas during the detector design process, in order to improve the performance 
for the measurements we described without reducing the effectiveness of the apparatus 
for other tasks. 
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